Concentration of solutions for a 
singularly perturbed neumann problem 
in non smooth domains 



Serena Dipierro 

SISSA, Sector of Mathematical Analysis 
Via Bonomea 265, 34136 Trieste, Italy 
E-mail address: dipierro@sissa.it 

Abstract 

We consider the equation — e 2 Ait + u — u p in a bounded domain Q, C R 3 with edges. 
We impose Neumann boundary conditions, assuming 1 < p < 5, and prove concentration of 
solutions at suitable points of dQ on the edges. 

1 Introduction 

In this paper we study the following singular perturbation problem with Neumann boundary 
condition in a bounded domain Q C M 3 whose boundary dQ is non smooth: 

— e 2 Au + u = u p in Q, 

§2=0 on dQ. ( > 

au 

Here p6 (1,5) is subcritical and v denotes the outer unit normal at dQ. 

Problem ([!]) or some of its variants arise in several physical and biological models. Consider, 
for example, the Nonlinear Schrodinger Equation 

dib fr 2 

lh Jt= ^ + V^- 1 \^- 2 ^ ) (2) 
at 2m 

where h is the Planck constant, V is the potential, and 7 and m are positive constants. Then 
standing waves of (|J) can be found setting ip (x, t) — e~ lEt / h v (x), where E is a constant and the 
real function v satisfies the elliptic equation 

- h 2 /\v + Vv= \v\ p ' 2 v 

for some modified potential V. In particular, when one considers the semiclassical limit h — > 0, 



the last equation becomes a singularly perturbed one; see for example | AM |, | FW |, and references 
therein. 

Concerning reaction-diffusion systems, this phenomenon is related to the so-called Turing's 
instability. More precisely, it is known that scalar reaction-diffusion equations in a convex domain 



admit only constant stable steady state solutions; see [CH], [Mat|. On the other hand, as noticed in 



1 



[ Tu I , reaction-diffusion systems with different diffusivities might generate non- homogeneous stable 
steady states. A well-known example is the Gierer-Meinhardt system, introduced in [GM| to 
describe some biological experiment. We refer to pj, [NTY| for more details. 



The study of the concentration phenomena at points for smooth domains is very rich and has 
been intensively developed in recent years. The search for such condensing solutions is essentially 
carried out by two methods. The first approach is variational and uses tools of the critical point 
theory or topological methods. A second way is to reduce the problem to a finite-dimensional one 
by means of Lyapunov-Schmidt reduction. 

The typical concentration behavior of solution Uq^ to ([!]) is via a scaling of the variables in 
the form 

U Q , e (x)~u(^), (3) 



where Q is some point of Q, and U is a solution of the problem 



AU + U = U P 



(or in 



{(X!,X 2 ,X 3 ) 



x 3 >0}), 



(4) 



the domain depending on whether Q lies in the interior of f2 or at the boundary; in the latter case 
Neumann conditions are imposed. When p < 5 (and indeed only if this inequality is satisfied), 
problem ([|) admits positive radial solutions which decay to zero at infinity; see |BL| ], ]&t[ . Solutions 
of (0) with this profile are called spike-layers, since they are highly concentrated near some point 
of fl 

Let us recall some known results. Boundary- spike layers are solutions of (Q) with a concentra- 
tion at one or more points of the boundary d£l as e — > 0. They are peaked near critical point of the 
mean curvature. It was shown in [NT1], [NT2 that mountain- pass solutions of ([!]) concentrate at 
dfl near global maxima of the mean curvature. One can see this fact considering the variational 
structure of ([!]). In fact, its solutions can be found as critical points of the following Euler-Lagrange 
functional 



~h{u) = - / (6 2 |V U 



dx 



\ p+1 dx, u e W 1 ' 2 (fi) . 



Plugging into I e a function of the form (||) with Q G dfl one sees that 

h (U Q , e ) = C e 3 - C^H (Q) + o (e 4 ) , 



(5) 



where Cq,C\ are positive constants depending only on the dimension and p, and H is the mean 
curvature; see for instance [AM], Lemma 9.7. To obtain this expansion one can use the radial 
symmetry of U and parametrize dfl as a normal graph near Q. From the above formula one can 
see that the bigger is the mean curvature the lower is the energy of this function: roughly speaking, 
boundary spike layers would tend to move along the gradient of H in order to minimize their energy. 
Moreover one can say that the energy of spike-layers is of order e 3 , which is proportional to the 
volume of their support, heuristically identified with a ball of radius e centered at the peak. There 
is an extensive literatur e regar ding t he sear ch of more genera l soluti o ns of (Q) concentrating at 
critical points of H; see f5FW|, @, jGPWj , [|iu|, S flLNTj , [|NPlfl , [jWel . 

There are other types of solutions of (|l|) with interior and/or boundary peaks, possible multiple, 
which ar e const r ucted by using gluing techniques or topological methods; see DW], [DY|, |GW |, 
| GW1 1 , [ G WW | , [Wa|. For interior spike solutions the distance function d from the boundary dil 
plays a role similar to that of the mean curvature H. In fact, solutions with interior peaks, as 
for the problem with the Dirichlet boundary condition, concentrate at critical points of d, in a 
generalized sense; see [LN|, [NW|, [Wei]. 
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Concerning a singularly perturbed problem with mixed Dirichlet and Neumann boundary con- 
ditions, in [ GMMPi] , ]GMMP2 | it was proved that, under suitable geometric conditions on the 



boundary of a smooth domain, there exist solutions which approach the intersection of the Neu- 
mann and the Dirichlet parts as the singular perturbation parameter tends to zero. 

There is an extensive literature regarding this type of problems, but only the case fl smooth 
was considered. Concerning the case non smooth, at our knowledge there is only a bifurcation 
result for the equation 

Am + Xf (u) = in ft, 

|^ = o on an, 



obtained by Shi in [3h| when 57 is a rectangle (0,a) x (0,6) in R 2 . 

In this paper we consider the problem (Q), where ft is a bounded domain in R 3 whose boundary 
dfl has smooth edges. If we denote by L an edge of dfl, we can consider the function a : T — > K 
which associates to every Q £ T the opening angle at Q, a (Q). As in the previous case, we can 
expect that the function a plays the same role as the mean curvature H for a smooth domain. In 
fact, plugging into I e a function of the form with Q £ T one obtains an expression similar to 
(|5|), with Coa (Q) instead of Co; see Lemma |4.3| . Roughly speaking, we can say that the energy of 
solutions is of order e 3 , which is proportional to the volume of their support, heuristically identified 
with a ball of radius e centered at the peak Q £ T; then, when we intersect this ball with the 
domain we obtain the dependence on the angle a (Q). 

The main result of this paper is the following 

Theorem 1.1. Let ft C R 3 be a piecewise smooth bounded domain whose boundary dft has a 
finite number of smooth edges, and 1 < p < 5. Fix an edge T, and suppose Q £ T is a local 
strict maximum or minimum of the function a, with a (Q) ^ n. Then for e > sufficiently small 
problem (Q) admits a solution concentrating at Q. 

Remark 1.2. The condition that Q is a local strict maximum or minimum of a can be replaced 
by the fact that there exists an open set V of T containing Q such that a (Q) > sup^y a or 
a (Q) < inigv «■ 

Remark 1.3. The condition a (Q) ^ n is natural since it is needed to ensure that dfl is not flat 
at Q. 

Remark 1.4. We expect a similar result to hold in higher dimension, with substantially the same 
proof. For simplicity we only treat the 3- dimensional case. 

The general strategy for proving Theorem |l.l| relies on a finite-dimensional reduction; see for 



example the book [AM] 

By the change of variables x ^ ex, problem (Q) can be transformed into 



—Am + u = u v in Q £ , 
= on <9fL 



§u _ n an (6) 

dv 

where Q e := iQ. Solutions of (@) can be found as critical points of the Euler-Lagrange functional 

^{U) = \j (|VM| 2 +M 2 )(fc l — [ \ U \ P+1 dx, M£V^ 1 ' 2 (a). (7) 

Now, first of all, one finds a manifold Z e of approximate solutions to the given problem, which 
are of the form Uq^ (x) — (p^ (ex) U (x — Q), where ip^ is a suitable cut-off function defined in a 



neighborhood of Q £ T; see the beginning of Section 4, Lemma 4.1 
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To apply the method described in Subsection 2.1 one needs the condition that the critical 
manifold Z e is non-degenerate, in the sense that it satisfies property ii) in Subsection 2.1. The 
result of non-degeneracy in r2 e , obtained in Lemma 4.2, follows from the non-degeneracy of a 



manifold Z of critical points of the unperturbed problem in K = K x K C R 3 , where K C K 2 
is a cone of opening angle a(Q). In fact, one sees that f2 £ tends to K as e — > 0. To show 
the non-degeneracy of the unperturbed manifold Z we follow the line of Lemma 4.1 in the book 
[A.M or Lemma 3.1 in [Ma|. We prove that A = is a simple eigenvalue of the linearized of the 
unperturbed problem at U S Z; see Lemma 3.1. Moreover, if a (Q) < n, it has only one negative 
simple eigenvalue; whereas, if a (Q) > n, it has two negative simple eigenvalues; see Corollary 



3.4. We note that in the case a (Q) — it, that is when dVL is flat at Q, A = is an eigenvalue of 
multiplicity 2. The proof relies on Fourier analysis, but in this case one needs spherical functions 
defined on a portion of the sphere instead of the whole S 2 . 

Then one solves the equation up to a vector parallel to the tangent plane of the manifold Z e , 
and generates a new manifold Z t close to Z e w hich represents a natural constraint for the Eulcr 
functional (|7|); see the proof of Proposition 4J3. By natural constraint we mean a set for which 
constrained critical points of I e are true critical points. 

We can finally apply the above mentioned perturbation method to reduce the problem to a finite 
dimensional one, and study the functional constrained on Z e . Lemma 4.3 provides an expansion of 
the energy of the approximate solution peaked at Q and allows us to see that the dominant term 
in the expression of the reduced functional at Q is a (Q). This implies Theorem 1.1. 



The paper is organized in the following way. In Section 2 we collect preliminary material: we 
recall the abstract variational perturbative scheme and obtain some useful geometric results. In 
Section 3 we prove the non-degeneracy of the critical manifold for the unperturbed problem in 
the cone K. In Section 4 we construct the manifold of approximate solutions, showing that it is 
a non-degenerate pseudo-critical manifold, expand the functional on the natural constraint and 
deduce Theorem 1.1. 



Notation 

Generic fixed constant will be denoted by C, and will be allowed to vary within a single line or 
formula. The symbols o e (1), or (1) o £j _r (1) will denote respectively a function depending on e 
that tends to as e -> 0, a function depending on R that tends to as R — ► +oo and a function 
depending on both e and R that tends to as e — > and R — > +oo. We will work in the space 
W 1 ' 2 (O e ), endowed with the norm ||w|| 2 = J Q (|Vu| 2 + u 2 ) dx, which we denote simply by ||it||, 
without any subscript. 



2 Some preliminaries 

In this section we introduce the abstract perturbation method which takes advantage of the vari- 
ational structure of the problem, and allows us to reduce it to a finite dimensional one. We refer 



the reader mainly to [AM], [Ma] and the bibliography therein. 

In the second part we make some computations concerning the parametrization of <9il and dfl e , 
and in particular of the edge. 

2.1 Perturbation in critical point theory 

In this subsection we recall some results about the existence of critical points for a class of func- 
tionals which are perturbative in nature. Given an Hilbert space H, which might depend on the 
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perturbation parameter e, let I e : H — > E be a functional of class C 2 which satisfies the following 
properties 

i) there exists a smooth finite-dimensional manifold, compact or not, Z t C H such that 
ll-^e( z )ll < Ce for every z £ Z t and for some fixed constant C, independent of z and e; 
moreover || I"(z) [q] \\ < Ce \\q\\ for every z£Z £ and every q £ T z Z e ; 

ii) letting P z : H — > (T z Z e ) , for every z G Z c , be the projection onto the orthogonal com- 
plement of T Z Z C , there exists C > 0, independent of z and e, such that P z I"(z), re- 
stricted to (T z Z e ) , is invertible from (T z Z e ) into itself, and the inverse operator satisfies 



(PzI'Jiz))- 1 



< c. 



We assume that Z e has a local C 2 parametric representation z = z%, £ £ R d . If we set W — (T z Z e ) , 
we look for critical points of I € in the form u = z + w with z £ Z e and w £ W. If P 2 : H — > W is 
as in m), the equation I' e (z + w) = is equivalent to the following system 

P z I' t (z + w) = (i/ie auxiliary equation) , , . 

(Id — P z ) I' e [z + w) = (the bifurcation equation) . 



Proposition 2.1. (See Proposition 2.2 in |Ma|/ j Let ZioZd. TTien t/iere exists eo > iwi/i i/ie 

following property: for all |e| < eo and for all z £ Z t , the auxiliary equation in (|J) /ias a unique 
solution w = w e (z) such that: 

j) w e (z) £ W is of class C 1 with respect to z £ Z e and w e (z) — > as |e| — > 0, uniformly with 
respect to z £ Z t , together with its derivative with respect to z, w' e ; 

jj) more precisely one has that \\w e (z)\\ = O (e) as e — »■ 0, for all z £ Z e . 

We shall now solve the bifurcation equation in (||). In order to do this, let us define the reduced 
functional <& e : Z e — > K by setting Q e (z) = I e (z + w e (z)). 



Theorem 2.2. (See Theorem 2.3 in jMa| /J Suppose we are in the situation of Proposition \2. l| ; and 
Zei ws assume that <!> e /ias, /or |e| sufficiently small, a critical point z e . Then u e — z e + w(z e ) is a 
critical point of / e ■ 



The next result is a useful criterion for applying Theorem [2J2| based on expanding I t on Z e in 
powers of e. 



Theorem 2.3. (See Theorem 2.4 in JMaj/ j Suppose the assumptions of Proposition 2.1 fto/d, and 
that for e small there is a local parametrization £ £ -U C R d o/ Z e such that, as e — > 0, I e admits 
the expansion I e (z^) = Co + eG(e£) + o(e), for £ £ -U, for some function G : U — > M. XTien 
we siiZZ /iave t/ie expansion ^> e (z^) = Co + eG(e£) -I- o(e), as e — > 0. Moreover, if £ £ U is a 
strict local maximum or minimum of G, then for |e| small the functional I e has a critical point u e . 
Furthermore, if £ is isolated, we can take u t — z^ e — o(l/e) as e — > 0. 

Remark 2.4. TTie last statement asserts that, once we scale back in e, the solution concentrates 
near £. 

2.2 Geometric preliminaries 

Let us describe <9f2 £ near a generic point Q on the edge T of dft e . Without loss of generality, 
we can assume that Q = £ M 3 , that xi-axis is the tangent line at Q to T in dil e , or dfl. In 
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a neighborhood of Q, let 7 : (— — >• M 2 be a local parametrization of T, that is (2:2,2:3) = 
7(2:1) = (71 (xi) ,72 (ail)). Then one has, for \xi\ < /i , 

(X2,X 3 ) = 7(2:1) 

= 7 (0)+7'(0) 2 ; 1 + i7"(0)x 2 + O(|x 1 | 3 ) 

= iy'(0)^+O(N 8 ). 

On the other hand, T is parametrized by (2:2,2:3) = 7 C (2:1) := ^7(6X1), for which the following 
expansions hold 



7e Ol) 



= | 7 "( )x 2 + O( e 2 |. 1 | 3 ), 
= e 7 "(0).Ti+O(e 2 |a;i| 2 ). 



(9) 



Now we introduce a new set of coordinates on Ban (Q) f] ^t- 

Vi=xi, (2/2,2/3) = (2:2,2:3) -j e (xi). 

The advantage of these coordinates is that the edge identifies with yi-axis, but the corresponding 
metric g ■ 
given by 

' 1 



metric g = (ff-y)^ will not be fiat anymore. If 7 £ (x\) = (7 e i (x\) , 7 e 2 (2:1)), the coefficients of g are 



(9ij) 



dx dx 



dyi dyi 

oi el 
dyi 

<)":,■ 2 



Of el 


dfe2 


dyi 


dyi 


1 








1 



From the estimates in (^) it follows that 

gii = Id + eA + 0(e 2 \ Xl \ 2 ) 

where 



(10) 



A = 





7" (Ofxi 



7" (0)2:i 




It is also easy to check that the inverse matrix (g' J ) is of the form g i: > = Id — eA + O (e 2 |xi| 2 ). 
Furthermore one has det g = 1. Therefore, by (SOI), for any smooth function u there holds 



A g u = Au — e 



(j" (0) VI ■ V te,,3)^J + (V (0) • ^(y 2 ,y 3) u) 

+0(e 2 |x 1 | 2 ) |V 2 U | + 0(e 2 | Sl | 2 ) |V«|. 



(11) 



Now, let us consider a smooth domain C M 3 and fl e — if2. In the same way we can describe 
d(l e near a generic point Q £ dfl e . Without loss of generality, we can assume that Q = £ R 3 , that 
{2:3 = 0} is the tangent plane of <9f2 e , or d(l, at Q, and that the outer normal v (Q) = (0,0, — 1). 
In a neighborhood of Q, let 23 = ip (xi,x 2 ) be a local parametrization of dQ,. Then one has, for 
| (2:1,2:2) I < Ml, 

2:3 = ip(xi,x 2 ) 

= \ ( A Q (xi,x 2 ) ■ (2:1,2:2)) + C Q (x 1 ,x 2 ) + 0(\ (2:1,2:2) | 4 ) , 



G 



where Aq is the Hessian of tp at (0, 0) and Cq is a cubic polynomial, which is given precisely by 

C Q (x 1 ,x 2 ) = - ^ dxidx dxk Mx iXjXk . 

i : j,k—l J 

On the other hand, dtl e is parametrized by x 3 = ip e (x\, x 2 ) := -i\> (exi, ex 2 ), for which the following 
expansions hold 

ip e (xi,x 2 ) = | (Aq (xi,x 2 ) • (xi,x 2 )) + e 2 Cq (x 1 ,x 2 ) + O (e 3 (xi, x 2 ) | 4 ) , 
^(x l7 x 2 ) = e(yiQ(x 1 ,x 2 )) j + e 2 i? Q (x 1 ,x 2 ) + 0( e 3 |(x 1 ,x 2 )| 3 ), (12) 
where Dq are quadratic forms in (xi,x 2 ) given by 

1 2 d 3 ib 

j,k=l J 



Concerning the outer normal v, we have also 



dip. 



' ~Bx^' ~ ^ ) ( i 
7 l = , , 9 = e(^Q (xi,x 2 )) +e 2 L» Q (xi,x 2 ) ,-1 + -e 2 |A Q (xi,x 2 ) 

+0( £ 3 |(x 1 ,x 2 )| 3 ). (13) 

Now we introduce a new set of coordinates on Bt± (Q) f] fa e - 

Zi=Xi, z 2 = x 2 , z 3 = x 3 - ip e (xi,x 2 ) . 

The advantage of these coordinates is that <9f2 e identifies with {23 = 0}, but, as before, the corre- 
sponding metric g = (gij)^ will not be fiat anymore. Its coefficients are given by 



\ 1 dV'e dip, dip, dip, dip, 

fir fir \ I ® z l dzi dzi dzo dz\ 

—— . I — dip^dip^ 2 J- dip^dip^ dip, 

~ \ f)v- Fiv ■ I ~~ ® z 2 Ozi dz 2 dz 2 dz 2 

dz\ dz2 



(fjij) 



From the estimates in ([12]) it follows that 

~g i3 ■ = Id + eA + e 2 B + 0(e 3 \( Zl ,z 2 )\ 3 ) , (14) 

where 



A = 

and 



A Q (zi,z 2 ) 
(A Q (z u z 2 )f 



p> _ ( A Q (^1^2) O Aq (Z1,Z 2 ) Dq(zi,Z 2 ) \ m 

\ (D Q (z u z 2 )) T 

It is also easy to check that the inverse matrix (<7 y ) is of the form g y = Id — eA + e 2 C 
0(e 3 |(zi,z 2 )| 3 ), where 



C = 



-D Q (z 1 ,z 2 ) 
{D Q (z 1 ,z 2 )) T \A Q {zi,z 2 )\ 2 



1 If the vector v has components (fi)^, the notation v ® v denotes the square matrix with entries (viVj) i y 
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Furthermore one has detg = 1. Therefore, by (^), for any smooth function u there holds 



AgU 



Au 



2 [D Q ■ V (Z1 



du 



\Aq (z l ,z 2 )\ 2 



du s 


) + trA Q 


du 
dz 3 


d 2 u 


— cUvDq 


du 


dz 3 dz 3 


dz 3 



+0 (e 3 | (z u z 2 ) | 3 ) |V 2 U | + O (e 3 | ( Zl ,z 2 ) | 3 ) |V«|. 



Moreover, from (13), we obtain the expression of the unit outer normal to dfl e , v, in the new 
coordinates z: 

/ 3 
e (A Q {z l7 z 2 )) + e 2 D Q (zi, z 2 ) , -1 + 2 £ Vq z *) \' 

+0(e 3 |(z 1 ,z 2 )| 3 ). 

Finally the area-element of dft e can be estimated as 

da = (l + O (e 2 (zi, z 2 ) | 2 )) dz\dz 2 . 

Now, locally, in a suitable neighborhood of Q G T, we can consider f2 as the intersection of 
two smooth domains f^i and Q 2 if the opening angle at Q is less than 7T, or as the union of them 

if the opening angle is greater than tt. In the first case one has dVl — (d&i D Cl 2 ^j U (j)Q 2 D , 
whereas in the second case dfl = D Cl 2 ^j U (d£l 2 D Q^j . Then, locally, one can straighten T 

and stretch the two parts of the boundary using the coordinates z for the smooth domains fii and 
2 - 



3 Study of the non degeneracy for the unperturbed problem 
in the cone 



Let us consider K = K x M. C R 3 , where K C R 2 is a cone of opening angle a, and the problem 

(15) 



Am + u = u p in K, 
on dK, 



g» = 



where p > 1 . 

If p < 5 and if u € I'K 1 ' 2 (K), solutions of ( |l5| ) can be found as critical points of the functional 
Ik : W 1 ' 2 (K) -> R defined as 



Zjf (u) = J / (|Vu| 2 + m 2 ) l — f \u\ p+1 dx. 

2 Jk P+ 1 iif 



(16) 



Note that 7^ is well defined on W 1 ' 2 (K); in fact, since K is Lipschitz, the Sobolev embeddings 
hold for p < 5; see for instance [Ad|, [ prifl . 

Let us consider also the elliptic equation in M 3 

-Au + u = u p , u G W 1 ' 2 (R 3 ) , u>0, (17) 

whi ch h as a positive radial solution U; see for instance AM |, [BL|, [Ma], ^tj. It has been shown 
in | Kw | that such a solution is unique. Moreover U and its radial derivatives decay to zero 
exponentially: more precisely satisfy the properties 



lim e r rU (r) = c 3 „, 

r— >+oo 



r U'(r) U"(r) 
hm = - hm = -1, 

i — s-+(xj (y (rj i — >+co U (r) 
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where r — \x\ and C3 iP is a positive constant depending only on the dimension n = 3 and p; see 



§3 



Now, if p is subcritical, the function U is also a solution of problem (|1 5[) - Moreover, if we 
consider a coordinate system with the si-axis coinciding with the edge of K, the problem (^) is 
invariant under a translation along the a?i-axis. This means that any 

U X1 (x) = U(x-(x 1 ,0,0)) 



is also a solution of (15). Then the functional I k has a non-compact critical manifold given by 

Z = {U X1 (x) : xi G R} ~ M. 

Now, to apply the results of the previous section, we have to characterize the spectrum and some 
eigenfunctions of I' K (U X1 ). More precisely we have to show the following 

Lemma 3.1. Suppose a G (0, 2ir) \ {tt}. Then the following properties are true: 

a) T Uxi Z = Ker [I« (U X1 )], for all x x G K; 

b) I'k (U Xl ) is an index Fredholm map ^ , for all X\ G M. 

Remark 3.2. The properties a) and b) imply that Z satisfies condition ii) in Subsection 2.1 and 
then it is non-degenerate for Ik- 

Proof. We will prove the lemma by taking x\ = 0, hence Uq = U . The case of a general x\ will 
follow immediately. 

Let us show a). It is known that there holds the inclusion TjjZ C Ker [I K (U)]; see for instance 
| AM), Section 2.2. Then it is sufficient to prove that Ker[I K (U)} C T V Z. Now, v G W 1 * 2 {K) 



belongs to Ker [I K (U)} if and only if 

-Av + v = pU p ~ l v in K 



|H = on QK. ^ 

We use the polar coordinates in K, r, 6, 95, where r>0,0<#<7r and < tp < a. Then we write 
v £ W 1:2 (K) in the form 

00 

v (xi,X2,x 3 ) = y^ v k (r) Y k (9, ip) , (19) 

where the Y k (9, (p) are the spherical functions satisfying 

-A S 2Y k = X k Y k mK, 

^ = p = 0,a. ( 2 °) 

Here A52 denotes the Laplace-Beltrami operator on S 2 (acting on the variables 6, <p). To determine 
Afc and the expression of Y k , let us split Y k as 



Y k (9, if) = J2 W (V) 



m=0 



2 A linear map T £ L (H, H) is Fredholm if the kernel is finite-dimensional and the image is closed and has finite 
codimension. The index of T is dim (Ker [T]) — codim (Ira [T]). 
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so that 



A S 2Y k = J2 
m=0 
oo r 

= E 



1 d ( . d_ 
smdde { de 



1 d 2 



'edip 2 



m=0 



1 



sin dd 



( S in06' fcim )$ 



k.m 



Then (po| ) becomes 

I K, m (o) = *U («) = o. 



1 Q (T)" 



Em=0^™^',™^^'" in if, 



(21) 



If we require that for all m 



-®'k m = Hm$k,m in [0, Of] , 



we obtain that <&k,m {f) = o-k,m cos {^ L{ p) satisfies (22) with fi r 
in (^lj) we have 



(22) 

Replacing this expression 



*U (0) = ^,™ («) = 0. 



-e fc . 



^fc,m ^^m— ^k,m®k,m^k,m in 



Since the &k,m are independent, we have to solve, for every m, the Sturm-Liouville equation 



1 



(sin 061 



k.m) 



sin 0d0 

Let us rewrite (^3|) in the following form 

1 d 



1 7r 2 m 2 



* k.m 



sin 



©fe,m — 0. 



(sin 06! 



1 7r 2 m z 



6fc,m — / ^fc6fc im , 



sin0 d0 N sin" ft a" 

so that we have to determine the eigenvalues Xk, m and the eigenfunctions of the operator 



(23) 



(24) 



1 



(sin 09' (0)) 



1 ir z m 



2^2 



■6(0). 



sin0G?0 v ~ - \ a ~ g - n 2 q a 2 
In order to do this, let us consider the case a = 7r, that is the following equation 



1 



sin d6 



i06' t 



1 2 

sir 



(25) 



Now, for every m, ( p5| ) has solution if Xk, m = k (fc + 1), with fc > \m\, and the solutions are the 
Legendre polynomials @k,m (0) = Pk,m (cos0); see for instance [Gro , [He], [Mu|, [ Mul |. Then, for 
a given value of k, there are 2k + 1 independent solutions of the form ®k,m (0) &k,m (</?), one for 
each integer m with — fc < m < k. Now, by the classical comparison principle, if we decrease a 
the corresponding eigenvalues Xk, m , given by (p4|), should increase, whereas if we increase a they 
should decrease; see for instance [Cha[. More precisely, if m = the equations ( p4[ ) and ( p5[ ) are 
the same, therefore the eigenvalues do not change (and they are 0,2,6,...). If m > 1 we cannot 
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give an explicit expression for the X k , m for general a, but we can use the comparison principle. In 
conclusion, we obtain that each Y k = J2m=o @fe,m^fc,m satisfies 

— A S 2Y k — Xk,mXk- (26) 

Now, one has that 

A (v k Y k ) = A r (v k ) Y k + ^v k A S 2Y k , (27) 

where A r denotes the Laplace operator in radial coordinates, that is A r — Jps + ^ -S^. Then, using 
©, (||) and (H), the condition (H) becomes 



E 



// 2 / Afc m TTV—1 

-v k v k + v k H —v k - pU p v k 



Y k = 0. 



fc=0 

Since the Y k are independent, we get the following equations for v k : 

2 \ 

Ak, m (v k ) := -v'l ~ -v' k + Vk + — jpVfc - pU p ~ x Vk = 0, m = 0, 1, 2..., k>m. 

Let us first consider the case m = 0. If k = 0, we have to find a vq such that 

2 

^o,o M = -v'o - -v' + v - pUP-^o = 0. 



It has been shown in [Kw[, Lemma 6, that all the solutions of ^4o,o ( v ) = are unbounded. Since 
we are looking for solutions vo G W 1,2 (M), it follows that Uq = 0. 
For k = 1 we have to solve 

2 2 
Ai.o (vi) = —v" - -v[ + vi + -rjVi - pU p ~ 1 vi = 0. 

Let U (r) denote the function such that U (x) = U (\x\), where U (x) is the solution of (p"7|). 
Reasoning as in the proof of Lemma 4.1 in [AM], we obtain that the family of solutions of Ai,o («i) = 
0, with vi £ W 1,2 (M), is given by V\ (r) — cU' (r), for some cel. 

Now, let us show that the equation Ak o ( u fc) = has only the trivial solution in W 1,2 (R), 
provided that k > 2. First of all, note that the operator A 10 has the solution U' which does not 
change sign in (0, oo) and therefore is a non-negative operator. In fact, if a denotes its smallest 
eigenvalue, any corresponding eigenfunction ip a does not change sign. If a < 0, then ip a should be 
orthogonal to U' and this is a contradiction. Thus a > and Ai$ is non- negative. Now, we can 
write 

A a _i_ A M ~ 2 

^fc,o — A i,o H ^2 • 

Since Afc.o — 2 > whenever fc > 2, it follows that A/^o is a positive operator. Thus A kt o (v k ) = 
implies that v k = 0. 

If m > 1 and a < n, using the comparison principle, we obtain that each X kjm is greater than 
2. Then, reasoning as above, we have that each v k = 0. 

Let us consider the case a > n. If m = 1 and fc = 1, using again the comparison principle, 
we have that < Ai.i < 2; whereas for m = 1, k > 2, and for m > 2, k > m, we have that 
each Afc.m > 2. Then in the last two cases we can use the non-negativity of the operator and 
conclude that v k — 0. In the case m = 1 and fc = 1 we note that the operator 

2 A 
Ai,i (vi) := -< - -ui + vi + -±-Vi - pU p - 1 v 1 
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has a negative eigenvalue, instead of the eigenvalue 0, since Ai.i < 2. Then also v\ = 0. 

Putting together all the previous information, we deduce that any v G Kcr [/" ({/)] has to be of 
the form 

v (xi,X2,x 3 ) = cU' (r) Y\ (0, ip) . 

Now, Yi is such that — A52 Y\ = X\mXx, namely it belongs to the kernel of the operator — A52 — 
Xi.rn.Id, and such a kernel is 1-dimensional. In conclusion, we find that 



v G span ^U'Yi^ = span /— — 1 = TjjZ. 



This proves that a) holds. It is also easy to check that the operator I K (U) is a compact perturbation 
of the identity, showing that b) holds true, too. This complete the proof of Lemma [O. □ 



Remark 3.3. Since U is a Mountain-Pass solution of (17), the spectrum of I K (U) has one 
negative simple eigenvalue, 1 — p, with eigenspace spanned by U itself. Moreover, we have shown 
in the preceding lemma that X = is an eigenvalue with multiplicity 1 and eigenspace spanned by 
If a < ir the rest of the spectrum is positive. Whereas if a > ir there is an other negative 

simple eigenvalue, corresponding to an eigenfunction U given by 



U (r, 9,ip)=u (r) cos (-<p) 6 (0) 



where satisfies ( |23[ ) with m = 1 and k = 1, and u satisfies the equation 

-v" - ^v' + v + ^-v- P U p - 1 v = 0. (28) 

From ( pq ) one has that there exists a positive constant C such that, for r sufficiently large, u (r) < 
Ce~* "I . In conclusion, one has the following result: 



Corollary 3.4. Let U and U be as above and consider the functional Ik given in (p_6[). Then for 
every x\ G R, U X1 (x) = U (x — {x\, 0, 0)) is a critical point of Ik- Moreover, the kernel of I K (U) 
is generated by . If a < tt the operator has only one negative eigenvalue, and therefore there 
exists S > such that 

r\j r 

I K {U)[v,v] > 8\\v\\ 2 , for every v G W 1 - 2 (K) , v±U, — . 

If a > it the operator has two negative eigenvalues, and therefore there exists 6 > such that 

dU 

VI (U) [v, v]>6 \\v\\ 2 , for every v G W 1 ' 2 (K) , v±U, U, — . 

ox 1 



4 Proof of Theorem 

For every Q on the edge T of <9f2 e , let fi = min {/Xj}, so that in Be. (Q) f] f2 e we can use the new 
set of coordinates z. Now we choose a cut-off function (p^ with the following properties 

(fin, (x) = 1 inB| (Q) , 

<Pfj, (x) =0 in R 3 \ Est (Q) , (29) 

|V^| + |V 2 ^I<C in B,(Qj\B,(Q). 
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For any Q £ F, we define the following function, in the coordinates (z\, z%, 23), 

U Q , t (z) := ^ (ez) U Q (z) , (30) 
where Uq (z) — U (z — Q). Then we consider the manifold 

z t = {u Q , e : q e r} . 

Now, we estimate the gradient of I e at Uq, £ , showing that Z e constitute a manifold of pseudo-critical 
points of I e . 

Lemma 4.1. There exists C > such that for e small there holds 

\\I'e(U Q ,e)\\ <Ce, forallQer. 



Proof. Let v £ W 12 (Cl e ). Since the function Uq, e is supported in B := Bj±_ (Q), see (p0[), we can 
use the coordinate z in this set, and we obtain 

I'e(U Q ,e)[v] = [ ^^vda+ f (-A- g U Q , e + U Q , e -\U Q , e \P)vdV s (z) 
= I + 11. 



Let us now estimate /: 

/=/ ^vda 1+ [ ^vda 2 =I 1+ I 2 . 
Jdsi el dvi J 9nc2 dv 2 

If K = K a (Q) denotes the cone of angle equal to the angle of the edge in Q, we have 

h = I (Uq (z) V<^ (ez) -vi + tp^ (ez) VUq (z) ■ v^vda-i 

JdK 

U Q (z) V</? m (ez) ■ (e (A Q (z l ,z 2 )) + e 2 D Q ( Zl ,z 2 ) , -1 + \e 2 \A Q ( Zl ,z 2 ) \ 2 

OK \ 1 

(ez) VU Q (z) ■ (e (A Q (z 1} z 2 )) + e 2 D Q ( Zl ,z 2 ) , -1 + \^\Aq *a) \ 2 

v(l + 0(e 2 \( Zl ,z 2 )\ 2 ))d Zl dz 2 
= a + b. 

Since \7(p^ (e-) is supported in K 3 \ Bj±_ (Q) and Uq has an exponential decay, we have that, for e 
small, 



\a\ < Cee~~ / \v\d Zl dz 2 . (31) 

JdK 



On the other hand 
b = 



■<\z-Q\<i 



-Q\<-1 



<p„ (ez) VU Q (z) ■ (e (A Q ( Zl ,z 2 )) + e 2 D Q (z u z 2 ) , -1 + \<?\Aq ( Zl ,z 2 ) \ 2 ^j 

v(l + 0(e 2 \(z 1 ,z 2 )\ 2 ))dz 1 dz 2 
<p M (ez) VUq (z) • ^ (A Q (z 1]Z2 )) + e 2 D Q ( Zl , z 2 ) , -1 + \e 2 \A Q ( Zl ,z 2 ) \ 2 ^j 

v(l + 0(e 2 \(y ll y 2 )\ 2 ))dy 1 dy 2 
<Cee~£ \v\dz x dz 2 + Ce / |W Q | • \v\dz x dz 2 . (32) 

JdK JdK 
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The estimates ( [3T| ) and (32), and the trace Sobolev inequalities imply |ii| < Ce ||u||. In the same 
way we can estimate I2, getting 

\I\<Ce\\v\\. 

Now let's evaluate II. Using (fTD) one has 



II 



K 



(-At/Q,e + U Q , e ~ \U Q , e \ P ) vdVg (z) 



K 



2 (7" (0) zi • V (22 , Z3) ^) + ( 7 " (0) • V {z2 , Za) U Q , e ) 



VdVg (z) 



+0(e 2 ) f (\ Zl \ 2 \V 2 U Q , e \ + \ Zl \ 2 \VU Qte \)vdV g (z) 



(33) 



= Ih + ell 2 + O (e 2 ) 7/3. 

Since AC/q iC = C/qA^ (ez) + 2VUq ■ V<^ (ez) + ^ (ez) A/7q and both Aip^ (e-) and V(^ M (e-) are 
supported in R 3 \ B±_ (Q), we get 



Hi = 



■<\x-Q\<i 



(-UqA^ (ez) - 2VU Q • V<^ (ez)) v (1 + O (e\z\)) dz 



fi<|z-QI<£ 



(-^ (ez) AC/q + U Q , e - \U Q A P ) v(1 + (e|z|)) dz 



\*-Q\<£ 



(-AU Q + U Q - \U Q \ P ) v(l+0 (e|z|)) dz. 



(34) 



Since Uq is a solution in R 3 the last term in ( |34| ) vanishes, and using the exponential decay of C/q 
at infinity and the properties of the cut-off function, see (p9|), one has 



\Ih \ < Ce~£ [ \v\dz. 
Jk 



By ( |30| ) we can compute also V( Z2>23 ) d g^' and ^ ( Z2iZ3 )Uq^ and we have 



Us 



2 7 " (0) zi 



K 



V 



( z 2,«3) " 



dip^iez) 
dzi 



dU, 



Q 



-2 7 " (0) Z! 



dip^ (ez) 



)^Q + ¥V (ez)V ( 

22,23.) 



dz\ 
dU Q 



-7" (0) • [V (Z2 , Z3) ^ (ez) U Q + <p„ (ez) V^^t/g] vdV g (z) 



<|z-QI<£ 



2 7 " (0) zi 



dtPfi (ez) 



g^Q . diff, (ez) , 



(z2 ' Za)_ dz^ — ^ Q + V ( Z2 ' Z3 )^ ^ + ~~dz~ 1 — V (^ 3 )Uq 



+7" (0) • V (Z2)Z3) ^ (ez) t/ Q ^ (z) 



2-QI<£ 



¥V ( ez ) 



27" (0) Zl ■ V(, 2 , Z3 )^f + 7" (0) • V (Z2 , Z3) C/ Q 



i>g?V^ (z) . 



Hence 



|/J 2 |<C /" 



<|z-QI<£ 



2|7"(n i |.i: l |(|C/ ( 3| + |^| + |V ( ^3 ) C/Q|)7h"(n,i. r,,| 



k-QI<j 



2|^ (ez) I -sup |t" (0)| N-|v ( 



dUQ, 

' dzi 1 



\v\dVg (z) 



+ |V (z2 , Z3 )^qI l«|d^(«) 
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Using again the exponential decay oHIq at infinity one can estimate the first term by Ce £ J K \v\dz 
and conclude that the second term is bounded. In the same way we can estimate J/3, getting 

|//|<CeH|. (35) 

From (133) and (B5l) we obtain the conclusion. □ 



Now, we need a result of non-degeneracy, which allows us to say that the operator I'J (Uq i6 ) is 
invertible on the orthogonal complement of Tjj q t Z IL . 

Lemma 4.2. There exists 6 > such that for e small, if a < ir, there holds 
l'l (U Q , e ) [v,v] > 5 \\vf , for every v € W 1 ' 2 (O e ) , v±U Q , £ , 
and, if a > 7T 7 there holds 

l' e '(U Q>e ) [v,v] >8\\vf, for every v G W 1 - 2 (a) , v±U Q , e , U Q ,e^f, 
where Uq^ is defined as Uq^ in (|30|). 

Proof. Let us consider the case a < ir. Let R ;§> 1; consider a radial smooth function xr : R 3 — > R 
such that 

Xr (x) = 1 in B R (0) , 

XR (x) =0 in R 3 \ B 2R (0) , (36) 
|V X r|<1 in B 2R (0) \ B R (0) , 

and set 

vi (x) = XR ( x ~ Q) v (x) , v 2 (x) = (1-Xr(x- Q)) v (x) . 
A straight computation yields 

IN 2 = IM| 2 + IM| 2 + 2 / (Vui • Vv 2 + viv 2 ) dx. 

We write j n (Vi>i • Vv 2 + viv 2 ) dx = 71 + 72, where 

71 = / Xr(1~Xr)(v 2 + \Vv\ 2 )dx, 

72 

Since the integrand in 72 is supported in B 2 r (Q) \ Br (Q), using (|36|) and the Young's inequality 
we obtain that I72I = or (1) ||u|| . As a consequence we have 

lkll 2 = lki|| 2 + ll« 2 || 2 + 27i + fl(i)lkll 2 . 

Now let us evaluate I'J {Uq^) [v, v] — o\ + a 2 + 03 , where 

°i = I" (U Q , e ) , o 2 = l'l (U Qje ) [v 2 ,v 2 ] , a 3 = 2I' e ' (U Q , e ) [v u v 2 ] . 



/ {y 2 Vv ■ \?xr - viVv ■ Vxr ~ v 2 \Vxr\ 2 ) dx. 
Jn, 
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Similarly to the previous estimates, since Uq decays exponentially away from Q, we get 

02 > C- 1 ^ 2 +o e , R (l)\\v\\ 2 , 

g 3 > C- X 7 i + o e , R (1) \\vf . (37) 
Hence it is sufficient to estimate the term o\ . From the exponential decay of Uq and the fact that 

dQ 



v±Uq^, dJ in' e it follows that 



( v ^ U Q,e) w i,2 {Qe) = - (V2, UQ, e ) w i,2(n c) = €<R (1) ||«|| , 

Vx,^A -U^A = ^(1)IM| 2 . (38) 



®Q /w 12 (o £ ) V ®Q / w 1 - 2 ^,) 



Moreover, since Uq^ is supported in B := Bj±_ (Q), see (|30|), we can use the coordinate z in this 
set, and we obtain 

( v i,UQ,e) w i,2 (n ) = [ v 1 dl i^' e vda+ [ vx(-AgUQ, e + U Qie )dVg(z) 
Jan, av Ja, 

= («i,^Q)w^ W +o £ (l)lhll, (39) 

where K — K a is the cone of opening angle equal to the angle of T in Q. In the same way we can 
obtain that 

/ ^ dJ b¥) = (*>1&) +*(i)IMI- ( 4 °) 



From the estimates ([38|), ([59|) and (ffOj), we deduce that for R sufficiently large and e sufficiently 
small 

( v uU Q ) wl , 2{K) = o e , R (l)\\vi\\, 

»i.§7f) = O e , R (l)\\ Vl \\. 

°W / W l ~ 2 (K) 

Now we can apply Lemma |3.1| getting 

I" ( u q) [vi,vi] > S \\vi\\ w i,2 {K) + o e , R (1) . 

Then the following estimate holds 

ax = I"(U Q )[v 1 ,v 1 ]+o e (l)\\v 1 \\>S\\v l \\ wl , 2(K) +o e , R (l)\\v\\ 

> <5|H| + 0e , fl (l)H. (41) 

In conclusion, from ( |37| ) and (|4l]) we deduce 

I' e '(U Q , £ )[v,v}>6\\v\\+o £ , R (l)\\v\\> 6 -\\v\\, 

provided R is taken large and e sufficiently small. This concludes the proof. 

The case a > n has substantially the same proof, but we h ave to consider also the function U 
and use the exponential decay of u at infinity, see Remark |3.3j . □ 

The following lemma provides an expansion of the functional I e (Uq^) with respect to Q. 
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Lemma 4.3. For e small the following expansion holds 

I e (U Q>e ) = Cva (Q) + 0(e), (42) 

where 

co= G-FTr)rr |tr,j(rrwjw9 - 



Proof. Since the function J7q j£ is supported in := Sji. (Q), see (30), we can use the coordinate 
z in this set, and we obtain 

I* (P Q , e ) = \( (\V g U Q J 2 + U* tC ) dV g (z)--±- I \U Q , e \ p+1 dV~ g (z) . 

Integrating by parts, we get 

Ie (u Q , e ) = - / u Q , e — + - / U Q , e (~A g U Q , e + U Q , e ) dV g (z) 
z JBndn c ° v ^ JBnn e 



1 



p + 1 



\u Q . t \ p+1 dv g {z) 

= 1 + 11, 



where / is the surface integral over the boundary and II refers to the last two terms. Now, / can 
be split in two terms which correspond to the surface integrals on the "faces" of the edge T: 



T 1 f JJ dU Q^~ M 1 f JJ dU Q'^~ ■ T 

1 = o / U Q>£ a- d(Tl + o / U Q>£ fl~ da2 = 7l + /2 - 

1 JBnan el av i z J Bnan t2 ° v i 



It is sufficient to evaluate I±, since the estimate of I2 is similar. Using the expression of /7q i£ , see 
©, we get 



h = ~ / U Q , e (C/qV^ (ez) + ^ (ez) VUq) 
1 -'Bnanji 



[e (A Q ( Zl ,z 2 )) + e 2 D Q (z 1: z 2 ) , -1 + ^e 2 \A Q (z u z 2 ) \ 2 ^j (l + O (e 2 \ (z u z 2 ) | 2 )) dz x dz 2 

= \\ <p„ {ez) C/qV^ m (ez) 

(e (Aq { Zl ,z 2 j) + e 2 D Q {zi,z 2 ) , -1 + \^Mq «a) | 2 ) (l + O (e 2 \ (z 1: z 2 ) \ 2 )) dz x dz 2 



1 



^ (ez) t/gVf/g 

k-QI<£ 



• (A Q (Z!,z 2 )) + e 2 D Q (z u z 2 ) , -1 + \^ 2 \Aq (zi,z 2 ) \ 2 ^j (l + O (e 2 \ (z Xl z 2 ) | 2 )) dz x dz 2 . 

Similarly to the previous estimates, we get I\=0 (e - ^) + O (e). Then we obtain that 

I = 0(e). (43) 
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Now, we have to evaluate 77: 



11= - 



+ - 



We have 



/ u Q , 

JBnn, 



Bnn c 



U Q , e (-AU QtC + U Q , e ) (1 + (e\z\)) dz 



2 7 " (0) Zl ■ V (Z2 , Z3) — 9± + 7 " (0) • V {Z2 , Z3) U Q 



II 



1 1 



2 p + 1 



+0(e 2 \ Zl f)- 



a(Q) 



1 



(l + 0(e\z\))dz 
\U Q , e f +1 (l + 0(e\z\))dz. 



P+ 1 J Etna 

rr 

\Uq (r) r +1 rsin 2 6drd6 + O (e) 



Putting together (|43|) and (44), we obtain ([42j) and this concludes the proof. 



(44) 
□ 



Let Pq : W 1 ' 2 (f2 e ) — > (Tu Q e Z < ^J ± be the projection onto the orthogonal complement of 
Tfj Q c Z e , for all Q on the edge T of dfl e . According to the lemma 4.2, we have that for e sufficiently 
small the operator Lq = Pq o I" (Uq i€ ) o Pq is invertible and there exists C > such that 



< C. 



Now, using the fact that I" (Uq^) is invertible on the orthogonal complement of Tjj q e Z t , we will 
solve the auxiliary equation. 

Proposition 4.4. Let I e be the functional defined in (R). Then for e > small there exists a 
unique w — w (e,Q) £ (Tuq e Z e ) such that I' e (L/q. £ + w (e, Q)) € T\j Q e Z £ . Moreover the function 
w (e, Q) is of class C 1 with respect to Q and there holds 



\\w(e,Q)\\<Ce, 



dw (e, Q) 



dQ 



< Ce. 



(45) 



Proof. We want to find a solution w € (Tu Q c Z^\ of PqI'^ [Uq^ + w) = 0. For every w G 
(Tu Q c Z e ) we can write 

I' e (U Q , t +W) = I' e (U Q , e ) + I'i (U Q , e ) [W] + R Q , e (W) , 

where Rq,e {w) is given by 

R Q , e (w) = I' e (U Q , t + w) - I[ (U Q . e ) - I'i (U Q , e ) [to] . 
Given v £ W 12 (fl e ) there holds 

RQ, e (w)[v] = - [ (\U Q , e + w\P-\U Q , e \ p -p\U Q , t \ p - 1 w)vdx. 



Using the following inequality 



| (a + b) p - a? - paP- l b\ < j 



C(p)\b\P 



for p < 2, 



C (p) (\b\ 2 + \b\P) forp>2, 
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for a, b 6 R, \a\ < 1, the Holder's inequality and the Sobolev embeddings we obtain 

\\Rq,< H Mil < C f (| W | 2 + \w\?) \v\dx < C (\\w\\ 2 + \\w\\ p ) \\v\\ . (46) 



Similarly, from the inequality 

\(a + b 1 f -(a + b 2 f -paP- 1 (b 1 -b 2 )\ 

/ C (p) (\h\P- 1 + \b2\P- 1 ) \h - b 2 \ forp<2, 
"\ C(p)(\b 1 \ + \b 2 \ + \b 1 \P- 1 + \b 2 \P- 1 )\b 1 ^b 2 \ forp>2, 

for a, bi, b 2 s K, |a| < 1, we get 

\\RQ,e (wi) [v] - Rq, c (w 2 ) H|| < C f (\wx\ + \w 2 \ + \wx\p- 1 + l^r 1 ) K - w 2 \ ■ \v\dx 

< c (\\ m \\ + ikii + n^ir 1 + n^ir 1 ) \\ Wl - W2 \\- \\v\\ . (47) 

Now, by the invertibility of the operator Lq — Pq o I" (Uq^) o Pq, we have that the function w 
solves PqI' € (Uq^ + w) = if and only if 

W = - (Lq)- 1 [P Q l' e (U Q , e ) + PqRq^ (w)\ . 

Setting 

Nq, £ (w) = - (LQ)' 1 [P Q l' e (U Q , e ) + PqRq, £ (w)] , 

we have to solve 

w = N Q ^ (w) . 

The norm of I' e (Uq^) has been estimated in Lemma Then from ([hj|) and ( |47| ) we obtain the 
two relations 

\\N Q , e (w)\\ < C ie + C 2 (|| W || 2 + |kf), (48) 

\\NQ, e {m)-NQ, e (wi)\\ < ^(iMl + lKll + lMr' + IKir 1 ) IK -w 2 || . (49) 

Now, for C > 0, we define the set 

W G = {w € {T UQ e Z e ) X : \\w\\<Ce}. 

We show that Nq^ is a contraction in Wq for C sufficiently large and for e small. Clearly, by ([48]), 
if C > 2Ci the set W c is mapped into itself if e is sufficiently small. Then, if w\, w 2 G W c , by 
©) there holds 

||iV Qie ( Wl ) - N Q>£ (w 2 )\\ < 2C (Ce + V- 1 ) \\ Wl -w 2 \\. 

Therefore, again if e is sufficiently small, the coefficient of \\wi — u>2|| in the last formula is less 
than 1. Hence the Contraction Mapping Theorem applies, yielding the existence of a solution w 
satisfying the condition 

\\w\\ < Ce. (50) 
This concludes the proof of the existence part. 
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Now the C 1 -dependence of the function w on Q follows from the Implicit Function Theorem; 
see also [AM], Proposition 8.7. In order to prove the second estimate in (45), let us consider the 
map H : iFx W 1 * 2 (O e ) xlxl — > W 1 ' 2 (Cl e ) x R defined by 



du Q ,. 



H (Q,w,a,e) = 



r e (UQ,e + w) - a-QQ 



Then w £ (T Uq € Z e ) is a solution of PqP £ (Uq 
£ W 12 (O e ) and /3 £ 

(Q,w,a,e)[v,f3] = 



Moreover, for v £ W 1 * 2 (O e ) and /3 £ R, there holds 
97? 



if and only if H (Q,w, a, e) 



d (w, a) 



r e '(u Q ,e + w) M-P^ 
r:(u Q , e )[v]-p2%f \ 



(51) 



+ o HI + IMI 



p-i 



To prove the last estimate it is sufficient to use the following inequality 



C(p) l&l"- 1 forp<2, 
C(p) (|6| + Ifep" 1 ) forp>2, 

for a,b £ R, |a| < 1, the Holder's inequality and the Sobolev embedding. Using the invertibility 
of the operator Lq = Pq o I" (£/g >e ) o Pq, it is easy to check that g^ 1 ^ {Q, 0, 0, e) is uniformly 
invertible in Q for e small. Hence, by ( |50|) and (0), also g^ 1 ^ (Q,w,a,e) is uniformly invertible 
in Q for e small. As a consequence, by the Implicit Function Theorem, the map Q <— > (wq,cxq) is 
of class C 1 . Now we are in position to provide the norm estimate of 9w qq^ ■ Differentiating the 
equation 

with respect to Q, we obtain 



H (Q,w Q ,a Q ,e) = 



OH 

ay o{w,a 



dH d(w Q ,a Q ) 
(Q,w,a,e) 



dQ 



Hence, by the uniform invertibility of g^ a ^ (Q, w, a, e) it follows that 



d(w Q ,a Q ) 



dQ 



< C 

< C 
C 



I>>(U Q , £ + w) 
(w. 

I' e '(U Q , £ + w) 
I'J (U Q;e + w) 



d 2 U Q ,, 
dQ 2 



u dQ 2 



dQ 
dU Q ,, 





+ IHI- 




) 


dQ 2 




dQ 2 





dQ 



+ \a\ + \\w\\ + e 



Note that a, similarly to w, satisfies \a\ < Ce. By the estimate in ( ^l| ) we obtain 



I'J (U Qtt + w) 



dQ 



< 



l"(U Q , e ) 



dU Q ,, 
dQ 



+ C[\\w\\ + \\w\\ 



p-i 
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Using the fact that 1" (Uq) 



9Uq 



I'J (Uq., + w) 



dQ 



< 



we obtain 
dUc 



For any v 6 W 1 ' 2 (K), one finds 

\(I' e '(U Q<e )-l"(U Q )) 



dzi 



, v 



I"(Uq) 



dUr 



Ce 



<pf \U Q<e -U Q \^-v + Ce. 



The last three formulas implies the estimate for ^"qq^ ■ This concludes the proof. 



□ 



Now we can state the following result, which allows us to perform a finite-dimensional reduction 
of problem (^) on the manifold Z e . 

Proposition 4.5. The functional >l/ c : Z e — > R defined by ^ e (Q) = I e (Uq^ + w (e, Q)) is of class 
C 1 in Q and satisfies 



K (Q) = o 



I' e (U Q ,e+w(e,Q)) = 0. 



Proof. This proposition can be proved using the arguments of Theorem 2.12 of AM]. From a 
geometric point of view, we consider the manifold 



Z e = {U Qie +w(e,Q):QeT}. 
Since (fill ) holds, we have that for e small 

Tu Q _ e Z e ~ T Uq e+w ( ey Q)Z e . 



(52) 



If Uq.c+w (e, Q) is a critical point of I e constrained on Z e , then I' e (Uq^ + w (e, Q)) is perpendicular 
to Tjj q e + w ( e ,Q)Z e , and hence, from (|52|), is almost perpendicular to Tu Q t Z t . Since, by construction 
of Z e , it is I' e (Uq^ + w (e, Q)) 6 T Uq e Z e , it must be I' e (Uq^ + w (e, Q)) = 0. This concludes the 
proof. 



□ 



4.1 Proof of Theorem [O 

First of all we have 

tt e (Q) = I e (U Q , e + w(e,Q)) 

= h (U Q<e ) + I' e (U Q>e ) [w (e, Q)) +0(\\w (e, Q)|f 



Now, using Lemma 4.1 and the estimate ([45|) we infer 

* e (Q) = I t (U Q , e ) + O (e 2 ) . 

Hence Lemma |4.3| yields 

* £ (Q) = C a (Q) + O (e) . 

Therefore, if Q € T is a local strict maximum or minimum of the function a, the thesis follows 
from Proposition |4.5| . 
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Abstract 

We consider the equation -e 2 Ati + u — u p in a bounded domain Q, C R 3 with edges. We 
impose Neumann boundary conditions, assuming 1 < p < 5, and we prove concentration of 
solutions at suitable points of dQ on the edges. 

1 Introduction 

In this paper we study the following singular perturbation problem with Neumann boundary 
condition in a bounded domain il C M 3 whose boundary dft is non smooth: 

|^ = o on an. 1 ' 

Here p G (1, 5) is subcritical and v denotes the outer unit normal at dQ. 

Problem ([!]) or some of its variants arise in several physical and biological models. Consider, 
for example, the Nonlinear Schrodinger Equation 

Bib h 2 

^^ = -7T^ + V^-lW- 2 ^, (2) 
at 2m 

where h is the Planck constant, V is the potential, and 7 and m are positive constants. Then 
standing waves of (^|) can be found setting ip (x, t) — e^ lEt l h v [x), where E is a constant and the 
real function v satisfies the elliptic equation 

- h 2 Av + Vv= \v\ p - 2 v 

for some modified potential V. In particular, when one considers the semiclassical limit h — > 0, 



the last equation becomes a singularly perturbed one (see for example | AM |, [EW], and references 
therein). 

Concerning reaction-diffusion systems, this phenomenon is related to the so-called Turing's 
instability. More precisely, it is known that scalar reaction-diffusion equations in a convex domain 



1 



admit only constant stable steady state solutions (see [OH], [Mat]). On the other hand, as noticed in 
[Tu|, reaction-diffusion systems with different diffusivities might generate non- homogeneous stable 
steady states. 

A well-known example is the Gierer-Meinhardt system, that is 



Ut = diAu 
v t — d 2 Av 

du 
dv - 



■ v H 
dv 

dv 



2y inOx(0,+oc), 
^ in Q x (0, +00) , 
= on dil x (0, +00) . 



(3) 



introduced in [GM] to describe some biological experiment. The functions u and v represent 
the densities of some chemical substances, the numbers p, q, r, s are non-negative and such that 
< < -xt, and it is assumed that the diffusivities d\ and d 2 satisfy d\ <C 1 <C d 2 . In the 
stationary case of (0), when d 2 — > +00 the function v is close to a constant (being nearly harmonic 
and with zero normal derivative at the boundary), and therefore the equation satisfied by u is 



similar to (pi), with e 2 = d%. We refer to (Ml, [NTY] for more details 



The study of the concentration phenomena at points for smooth domains is very rich and has 
been intensively developed in recent years. The search for such condensing solutions is essentially 
carried out by two methods. The first approach is variational and uses tools of the critical point 
theory or topological methods. A second way is to reduce the problem to a finite-dimensional one 
by means of Lyapunov-Schmidt reduction. 

The typical concentration behavior of solution u t to ([j]) is via a scaling of the variables in the 
form 

' x - Q~~ 



U Q ,e (a 



U 



where Q is some point of Q, and U is a solution of the problem 



- AU + U = U p 



(or in Rl = {(ssi.xa.a*) G K 3 : ^3 > 0}), 



(4) 



(5) 



the domain depending on whether Q lies in the interior of f2 or at the boundary; in the latter case 
Neumann conditions are imposed. When p < 5 (and indeed only if this inequality is satisfied), 
problem (g) admits positive radial solutions which decay to zero at infinity, sec [ BL ] , |St| . Solutions 
of (|l]) with this profile are called spike-layers, since they are highly concentrated near some point 
of Q. 

Let us recall some known results. Boundary-spike layers are solutions of (Q) with a concentra- 
tion at one or more points of the b ound ar y d£l as e — > 0. They are peaked near critical point of the 
mean curvature. It was shown in [ NT1 |, [ NT2 that mountain-pass solutions of ([!]) concentrate at 
dfl near global maxima of the mean curvature. One can see this fact considering the variational 
structure of (|j). In fact, its solutions can be found as critical points of the following Euler-Lagrange 
functional 



/.(«) = 5 



dx — 



p+1 



\ p+1 dx, u s w 1 ' 2 (n) . 



Plugging into I e a function of the form (^) with Q £ <9f2 one sees that 

h (u Q , e ) = C e 3 - de A H (Q) + o (e 4 ) , 



(6) 



where Cq,Ci are positive constants depending only on the dimension and p, and H is the mean 
curvature (see for instance AM ], Lemma 9.7). To obtain this expansion one can use the radial 
symmetry of U and parametrize d£l as a normal graph near Q. From the above formula one can see 



2 



that the bigger is the mean curvature the lower is the energy of this function: roughly speaking, 
boundary spike layers would tend to move along the gradient of H in order to minimize their 
energy. Moreover one can say that the energy of spike-layers is of order e 3 , which is proportional 
to the volume of their support, heuristically identified with a ball of radius e centered at the peak. 
There is an extensive literature regarding the search of more general solutions of (Q) concentrating 
at critical points of if; see pFW|] , pr| ]GPW[ , @, JCnt]], flNPTj ], jWeJ. 

There are other types of solutions of (111) with interior and/or boundary peaks, possible multiple, 
which arc constructed by using gluing techniques or topological methods; see |DW], [DY|, ]GW 



[ GW1 ] , [ GWW | , [ Wa | . For interior spike solutions the distance function d from the boundary dfl 
plays a role similar to that of the mean curvature H. In fact, solutions with interior peaks, as 
for the problem with the D irich lct bou ndary condition, concentrate at critical points of d, in a 
generalized sense; see [LN|, | NW |. 



Wei 



Concern ing a sing ul arly pertu rbed problem with mixed Dirichlet and Neumann boundary con- 
ditions, in [ GMMPij , ||GMMP2 | it was proved that, under suitable geometric conditions on the 
boundary of a smooth domain, there exist solutions which approach the intersection of the Neu- 
mann and the Dirichlet parts as the singular perturbation parameter tends to zero. 

There is an extensive literature regarding this type of problems, but only the case smooth 
was considered. Concerning the case fl non-smooth, at our knowledge there is only a bifurcation 
result for the equation 

Au + Xf (u) = in O, 
fi = on dn, 

au 



obtained by Shi in [3h| when f2 is a rectangle (0,a) x (0,6) in R 2 . 

In this paper we consider the problem (Q), where Q is a bounded domain in R 3 whose boundary 
dil has smooth edges. If we denote by T an edge of dtt, we can consider the function a : T — > M. 
which associates to every Q G T the opening angle at Q, a(Q). As in the previous case, we can 
expect that the function a plays the same role as the mean curvature H for a smooth domain. In 
fact, plugging into i c a function of the form ([I) with Q E T one obtains an expression similar to 
(H), with Coot (Q) instead of Co (see Lemma L3). Roughly speaking, we can say that the energy of 
solutions is of order e 3 , which is proportional to the volume of their support, heuristically identified 
with a ball of radius e centered at the peak Q G T; then, when we intersect this ball with the 
domain we obtain the dependence on the angle a (Q). 

The main result of this paper is the following 

Theorem 1.1. Let C R 3 be a piecewise smooth bounded domain whose boundary dfl has a 
finite number of smooth edges, and 1 < p < 5. Fix an edge T, and suppose Q € T is a local 
strict maximum or minimum of the function a, with a (Q) =/= n. Then for e > sufficiently small 
problem ([!]) admits a solution concentrating at Q. 

Remark 1.2. The condition that Q is a local strict maximum or minimum of a can be replaced 
by the fact that there exists an open set V of T containing Q such that a (Q) > sup^y a or 
a (Q) < inigv a. 

Remark 1.3. The condition a(Q) ^ n is natural since it is needed to ensure that dfl is not flat 
at Q. 

Remark 1.4. We expect a similar result to hold in higher dimension, with substantially the same 
proof. For simplicity we only treat the 3- dimensional case. 



The general strategy for proving Theorem 1.1 relies on a finite-dimensional reduction, see for 
example the book [AM]. 
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By the change of variables ii-yei, problem ([!]) can be transformed into 



-Au + u — u p in fi e , 
|s = on <9ft e 



(7) 



where f2 c = if2. 

Solutions of (0) can be found as critical points of the Euler-Lagrange functional 
I e (u) = l f (\Vu\ 2 + u 2 )dx l — [ \u\ p+1 dx, ueW 1 ' 2 ^). 



(8) 



Now, first of all, one finds a manifold Z t of approximate solutions to the given problem, which 
are of the form Uq^ (x) = (p^ (ex) U (x — Q), where ip^ is a suitable cut-off function defined in a 
neighborhood of Q G T (see the beginning of Section 4, Lemma 4.1). 

To apply the method described in Subsection 2.1 one needs the condition that the critical 
manifold Z e is non-degenerate, see (||) and (JTJJ). The result of non-degeneracy in f2 e , obtained in 
Lemma [T^, follows from the non-degeneracy of a manifold Z of critical points of the unperturbed 
problem in K — K x R C R 3 , where K C R 2 is a cone of opening angle a (Q). In fact, one sees that 
f2 € tends to K as e — »■ 0. To show the non-degeneracy of the unperturbed manifold Z we follow 
the line of Lemma 4.1 in the book [AM| or Lemma 3.1 in [Ma]. We prove that A = is a simple 



eigenvalue of the linearized of the unperturbed problem at U £ Z, see Lemma 3T. Moreover, if 



a (Q) < 7r, it has only one negative simple eigenvalue; whereas, if a (Q) > 7r, it has two negative 
simple eigenvalues. We note that in the case a (Q) = ir, that is when d£l is flat at Q, the A = 
is an eigenvalue of multiplicity 2. The proof relies on Fourier analysis, but in this case one needs 
spherical functions defined on a portion of the sphere instead of the whole S 2 . 

Then one solves the equation up to a vector parallel to the tangent plane of the manifold Z e , 
and generates a new manifold Z e close to Z e which represents a natural constraint for the Eulcr 
functional see the proof of Proposition 4.5. By natural constraint we mean a set for which 
constrained critical points of I t are true critical points. 

We can finally apply the above mentioned perturbation method to reduce the problem to a finite 
dimensional one, and study the functional constrained on Z t . Lemma 4.S provides an expansion of 
the energy of the approximate solution peaked at Q and allows us to see that the dominant term 
in the expression of the reduced functional at Q is a (Q). This implies Theorem 1.1. 



The paper is organized in the following way. In Section 2 we collect some preliminary material: 
we recall the abstract variational perturbative scheme and obtain some geometric results. In 
Section 3 we prove the non-degeneracy of the critical manifold for the unperturbed problem in 
the cone K . In Section 4 we construct the manifold of approximate solutions, showing that it is 
a non-degenerate pseudo-critical manifold, expand the functional on the natural constraint and 



deduce Theorem 1.1 



Notation 

Generic fixed constant will be denoted by C, and will be allowed to vary within a single line or 
formula. The symbols o e (1), 0^(1) o e ^n (1) will denote respectively a function depending on e 
that tends to as e — >• 0, a function depending on R that tends to as R — >• +oo and a function 
depending on both e and R that tends to as e — > and R —> +oo. We will work in the space 
W 1 ' 2 (fi £ ), endowed with the norm ||u|| 2 = J Q (|Vu| 2 + u 2 ) dx, which we denote simply by ||«||, 
without any subscript. 
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2 Some preliminaries 



In this section we introduce the abstract perturbation method which takes advantage of the vari- 
ational structure of the problem, and allows us to reduce it to a finite dimensional one. We refer 



the reader mainly to [AM], [Ma] and the bibliography therein. 

In the second part we make some computations concerning the parametrization of dfl (and 
dfl e ), and in particular of the edge. 

2.1 Perturbation in critical point theory 

In this subsection we will discuss the existence of critical points for a class of functionals which 
are perturbative in nature. Given an Hilbert space H (which might depend on the perturbation 
parameter e), let J e : H — >• K be a functional of class C 2 which satisfies the following properties 

i) there exists a smooth finite-dimensional manifold (compact or not) Z t C H such that 
ll-^eWII < Ce for every z £ Z e and for some fixed constant C (independent of z and e); 
moreover ||/"(z) [q]\\ < Ce \\q\\ for every z £ Z t and every q £ T z Z e ; 

ii) letting P z : H — > (T z Z e ) , for every z £ Z e , be the projection onto the orthogonal com- 
plement of T z Z e , there exists C > (independent of z and e) such that P z I"(z), re- 
stricted to (T z Z e ) , is invertible from (T z Z e ) into itself, and the inverse operator satisfies 



< c. 

For example, one can consider functionals of the form 

I e (u)=I (u) + eG(u), 

where Iq £ C 2,a (H, R) is the unperturbed functional and G £ C 2,a (H, K) is a perturbation. 
Suppose there exists a smooth <i-dimensional critical manifold of Iq, that is a manifold Z such that 
all z £ Z is a critical point of Iq. Let T Z Z denote the tangent space to Z at z. Since Z is a critical 
manifold then for every z £ Z one has Io( z ) = 0- Differentiating this identity, we get 

(I^(z)[v}\4>)=0, VveT z Z,Vcj)eH, 

and this implies that every v £ T Z Z is a solution of the linearized equation Iq(z) [v] — 0, namely 
that v £ Ker [Iq(z)]. Then one has T Z Z C Ker [Iq (z)]. In particular, ig (z) has a non-trivial kernel 
(whose dimension is at least d) and hence all the z £ Z are degenerate critical points of Iq. If we 
want that this degeneracy is minimal, we have to require that 

T Z Z = Ker [7 '(z)] , Vz £ Z. (9) 

In addition to (^|) we assume also that 

for all z £ Z, Iq (z) is an index Fredholm mapj^] (10) 

The last two properties imply that Iq(z) is invertible from (T z Z e ) onto itself. Therefore, if I(f 
and G are bounded (near Z) and if we take Z t = Z, properties i), ii) above will follow. 

Now, we assume that Z^ has a (local) C 2 parametric representation z = zt, £ £ M. d . If we set 
W — (T z Z e ) , we look for critical points of I e in the form u — z + w with z £ Z e and w £ W. If 



X A linear map T £ L (H, H) is Fredholm if the kernel is finite-dimensional and the image is closed and has finite 
codimension. The index of T is dim (Ker [T]) — codim (Im [T]). 
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P z : H — > W denotes the orthogonal projection onto W, the equation I' e (z + w) = is equivalent 
to the following system 

P z I' e (z + w) = (the auxiliary equation) , , . 

(Id — P z ) I' e (z + w) =0 (the bifurcation equation) . 



Proposition 2.1. (See Proposition 2.2 in fMa|/ j Let hold. Then there exists eo > with the 

following property: for all |e| < eo and for all z G Z t , the auxiliary equation in (11) has a unique 
solution w = w e (z) such that: 

j) W e (z) G W is of class C 1 with respect to z G Z e and w e (z) — > as |e| — > ; uniformly with 
respect to z £ Z e , together with its derivative with respect to z, w' e ; 

jj) more precisely one has that ||w £ (z)|| —0(e) as e -> 0, for all z G Z e . 



We shall now solve the bifurcation equation in (|1 1|) . In order to do this, let us define the reduced 
functional <I> C : Z e — >• R by setting 

$ £ (z) =/ e (z + w e (z)). 

From a geometric point of view the argument can be outlined as follows. Consider the manifold 
Z e — {z + W e (z) : z G Z e }. If z e is a critical point of $ e , it follows that u e = z e + w(z e ) G Z e is a 
critical point of J e con stra ined on Z e and thus u e satisfies I' e (u e ) _L T Ut Z t . Moreover, the definition 
of w e , see Proposition 2J, implies that J £ (z + w e (z)) G T 2 Z £ . I n pa rticular, I' e (u e ) G T Ze Z e . Since, 
for |e| small, T Uc Z t and T Zt Z e are close, see j) in Proposition |2.l| , which is a consequence of the 
smallness of w' t , it follows that I' e (u e ) — 0. A manifold with these properties is called a natural 
constraint for L. 



2.1 



and 



Theorem 2.2. (See Theorem 2.3 in t MaJ) Suppose we are in the situation of Proposition 
let us assume that <i> e has, for |e| sufficiently small, a critical point z e . Then u e = z e + w(z e ) is a 
critical point of I t . 



The next result is a criterion for applying Theorem 2.2, based on expanding I t on Z e in powers 
of e. 



Theorem 2.3. (See Theorem 2.4 in \Mq] ) Suppose the assumptions of Proposition 2.1 hold, and 
that for e small there is a local parametrization £ G -U C R d o/ Z e smc/i £/iat 7 as e — >• 0, I e admits 
the expansion 

I e (zs) = C + eG(eO + o(e), £ G -^U 

for some function G : U — > M. TTien we siiZi /iawe t/ie expansion 

®e(z() = C + eG(e£) + o(e), as e -)• 0. 

Moreover, if ^ G U is a strict local maximum or minimum of G, then for |e| small the functional 
I e has a critical point u e . Furthermore, if £ is isolated, we can take u e — z^ e = o(l/e) as e — > 0. 

Remark 2.4. The last statement asserts that, once we scale back in e, the solution concentrates 
near £. 



G 



2.2 Geometric preliminaries 

Let us describe Sf2 c near a generic point Q on the edge V of dft e . Without loss of generality, we 
can assume that Q = 6 I 3 , that Xi-axis is the tangent line at Q to T in <9f2 e (or <9fi). In a 
neighborhood of Q, let 7 : (— /io,/io) -4 I 2 be a local parametrization of T, that is (x2,xs) = 
7 (xi) = (71 (xi) ,72 (xi)). Then one has, for |xj.| < /Lt , 

(a; 2 ,x 3 ) = 7(»i) 

= 7 (0)+7'(0)x 1 + i7"(0)^ + O(|x 1 | 3 ) 

= \l" (0)4+0 {\ Xl \ a ). 

On the other hand, T is parametrized by (0:2, £3) = 7e (^i) := -j(exi), for which the following 
expansion holds 

7e(xi) = | 7 "(0)^ + O(e 2 | : r 1 | a ) 

^ = e7 "(0)x 1 +O( £ 2 | a ; 1 | 2 ). (12) 

Now we introduce a new set of coordinates on Bto. (Q) f] 

yi=xi; (j/2, 2/3) = (>2,Z3) -7e (zi) • 

The advantage of these coordinates is that the edge identifies with yi-axis, but the corresponding 
metric g — (gij)^ will not be fiat anymore. If 7e (xi) = ( 7e i (xi) , 7 e 2 (xi)), the coefficients of g are 
given by 







dx dx 



\ _|_ 87el <?7tl _|_ t?7e2 C>7e: 



d~f t l 


C?7e2 


dyi 


dyi 


1 








1 



From the estimates in ([12]) it follows that 

9ij = Id + eA + 0(e 2 \ Xl \ 2 ) , (13) 

where 

U"(o) T ^i )• 

It is also easy to check that the inverse matrix (<?' J ) is of the form g lJ — Id — eA + O (e 2 |xi| 2 ). 
Furthermore one has det g = 1. Therefore, by (|13|), for any smooth function u there holds 



A g u = Au — e 



2 (7" (0) yi • V — J + (7" (0) • V {y2 , ya) u) 

+0(6 2 | a;i | 2 )|V 2 U | + 0(e 2 |a ;i | 2 )|V M |. (14) 



Now, let us consider a smooth domain O g R 3 and f2 e = ifi. In the same way we can 
describe the <90 e near a generic point Q £ d£l e . Without loss of generality, we can assume that 
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Q = 6 R 3 , that {x3 = 0} is the tangent plane of dCl e (or dfl) at Q, and that the outer normal 
v (Q) = (0) 0, —1). In a neighborhood of Q, let X3 = tp (xi, x 2 ) be a local parametrization of c?f2. 
Then one has, for | (xi, X2) \ < /xi, 

x 3 = -0(xi, x 2 ) 

= 2 ( Xl ' X2 ) ■ (^I'^a)) + Cq Oi,^) +0 (I (Xi,x 2 ) | 4 ) , 

where ^4q is the Hessian of ^ at (0, 0) and Cq is a cubic polynomial, which is given precisely by 

C Q (x 1 ,x 2 ) = - dxid 9xk Mx iXj x k . 

i,j,k—l J 

On the other hand, dtl e is parametrized by £3 = ip e (xi, x 2 ) := ^V 1 ( e:c ii ea; 2), for which the following 
expansion holds 

i/> e (xi,x 2 ) = I (A Q (xi,x 2 ) • (xi,x 2 )) + e 2 C Q (xi,x 2 ) + O (e 3 (xi,x 2 ) | 4 ) , 
|^(xi,x 2 ) = e (^ Q (x 1 ,x 2 )) i + e 2 ^(xi,x 2 ) + 0( e 3 |(x 1 ,x 2 )| 3 ), (15) 

where Dq are quadratic forms in (xi,x 2 ) given by 

2 

2 ^ (>.r,().r ,(lv] 

Concerning the outer normal ^, we have also 

" = 7= ^t ^tT = U(AQ(x 1 ,x 2 ))+e 2 C Q ( a ;i,a; 2 ),-l + -e 2 | J 4Q(xi,x 2 )| 2 

+0( £ 3 |(x 1 ,x 2 )| 3 ). (16) 

Now we introduce a new set of coordinates on Bn (Q) f] f) e . 

z 1 =x 1 ; z 2 = x 2 ; z 3 = x 3 -ip e (x 1 ,x 2 ). 

The advantage of these coordinates is that dCl e identifies with {z 3 } = 0, but the corresponding 
metric g = {fjij)^ will not be flat anymore. Its coefficients are given by 



1 ®' i '4 ) 
D i Q {xi,x 2 ) = - J2 aww^ (0,0)gjJfc- 



/ /fl-r r)T \ / dzi dzi dzi dz2 dz\ 

(~ \ — I _ __ I — djKdjK 1 1 JiVk M<l 

WJV - I q ' a \ dz x d Zl 1 + dz 2 dz 2 dz 2 

\ UZ l UZ J / \ g»/> £ 1 



021 dz2 



From the estimates in (|l5]) it follows that 

5y =Id + eA + e 2 B + 0(e 3 \{z llZ2 )\ 3 ), (17) 

where 



.4 



A Q (z 1 ,z 2 ) 



(A Q (z 1 ,z 2 )f 



8 



and 



B 



Aq (zi,z 2 ) ® A Q (zi,z 2 ) D q (zi,z 2 ) 
{D Q (z u z 2 )) T 

It is also easy to check that the inverse matrix ) is of the form g l i 
(e 3 | (zi,z 2 ) | 3 ), where 



Id - eA + e 2 C + 



C 



-D Q {zi,z 2 ) 



-{D Q {z x ,z 2 )Y \A Q ( Zl ,z 2 )\ 2 
Furthermore one has detg = 1. Therefore, by (0), for any smooth function u there holds 



AgU = Am — e 



-2 D. 



„ / . , . _ du \ du 
2 \Aq (z u z 2 )-V (zuZ2) —\ + trA Q — 

du\ 2 <9 2 it 

■ |-4q (21^2) | 2 



divD 



du 
dz 3 

+0 (e 3 \ (z u z 2 ) | 3 ) |V 2 U | + O (e 3 | z 2 ) | 3 ) |V«|. 



Moreover, from (16), we obtain the expression of the unit outer normal to dfl e , z/, in the new 
coordinates z: 

^ = ( e ( A Q ( Zl > Z2 )) + e2 °Q ( 2l > Z2 ) ' ~ X + l e2 \ A Q ( z ^ 7 ^) \' 

+0(e 3 |(z 1 ,z 2 )| 3 ). 

Finally the area-element of dCt e can be estimated as 

da = (l + O (e 2 (zi, z 2 ) | 2 )) dz\dz 2 . 

Now, locally, in a suitable neighborhood of Q s T, we can consider as the intersection of 
two smooth domains f^i and Sl 2 if the opening angle at Q is less than tt, or as the union of them 
if the opening angle is greater than it. In the first case one has dfl = (d&i n Cl^ U (dfl 2 n Clij , 

whereas in the second case (d£l 2 (iCl^j. Then, locally, one can straighten T 

and stretch the two parts of the boundary using the coordinates z for the smooth domains fii and 

n 2 . 



3 Study of the non degeneracy for the unperturbed problem 
in the cone 

Let us consider K = K x R C R 3 , where K C R 2 is a cone of opening angle a, and the problem 



-Am + u = u p in K, 
= on dK, 



(18) 



where p > 1. 

If p < 5 and if u S VF 1,2 (if), solutions of (|l8| ) can be found as critical points of the functional 
I K : W 1 ' 2 (K) -> R defined as 



J* («) = J / (|Vm| 2 + M 2 ) ' 



p+ 1 



| p+1 <iE. 



(19) 



A 



2 If the vector u has components (t>;)^, the notation i> <%> v denotes the square matrix with entries (viVj) . 
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Note that Ik is well defined on W 1 ' 2 (K); in fact, since K is Lipschitz, the Sobolev embeddings 
hold for p < 5 (see for instance [Ad], [Gri]). 

Let us consider also the elliptic equation in M 3 



Au + u — u p . 



G W 1 ' 2 (R 3 ) , u > 0, 



(20) 



which has a positive radial solution U (see for instance |AM |, [BL], [Ma], pt|). It has been 
shown in [Kw] that such a solution is unique. Moreover U and its radial derivative decay to zero 
exponentially: more precisely satisfy the properties (see [BL]) 

U'(r) . U"(r) 
lim e r rlf (r) = C3 „; lim — — - = — lim — - = —1, 

where r = \x\ and C3. p is a positive constant depending only on the dimension n = 3 and p. 

Now, if p is subcritical, the function U is also a solution of problem (|l8|). Moreover, if we 
consider a coordinate system with the xi-axis coinciding with the edge of K, the problem (|l^) is 
invariant under a translation along the Xi-axis. This means that any 

U X1 (x) = U(x-( Xl ,0,0)) 



is also a solution of (18). Then the functional Ik has a non-compact critical manifold given by 

Z = {U X1 (x) : x x e K} ~ M. 

Now, to apply the results of the previous section, we have to characterize the spectrum and some 
eigenfunctions of 1'^ (U Xl )- More precisely we have to show the following 

Lemma 3.1. Suppose a G (0,27r) \ {it}. Then Z is non-degenerate for Ik, namely the following 
properties are true: 

i) T Uxi Z = Ker [I' K (U Xl )\, Vsi G R; 

ii) I' K (U Xl ) is an index Fredholm map, for all X\ G R. 



Proof. We will prove the lemma by taking x\ = 0, hence C/o = U . The case of a general x\ will 
follow immediately. 

Let us show i). We have seen in Subsection 2.1 that there holds the inclusion T\jZ C 
Ker [I K (U)]. Then it is sufficient to prove that Ker [I K (U)] C T V Z. Now, v G IV 1 ' 2 {K) be- 
longs to Ker [I' K ([/)] if and only if 

-Av + v — pU p ~ l v in K, , . 

fi = on dK. ^ ' 

We use the polar coordinates in K, r, 6, ip, where r>0, < 6 < tt and < <p < a. Then we write 
v G W 1 - 2 (K) in the form 

oo 

V(X!,X2,X 3 ) 

where the Yk (0, tp) are the spherical functions satisfying 

f -A S 2Yk = X k Yk in K 



10 



Here A52 denotes the Laplacc-Bcltrami operator on S 2 (acting on the variables 9, f). To determine 
Afc and the expression of Y k , let us split Y k as 



Y k {9,if) = ®k,m(0)$k,m (V) 



so that 



E 

m— 
00 r 

E 

m=0 



m=0 



1 3 ( . n d 
sm 



sin 6*961 V 56* ) sin 2 9 dip 2 



1 a 2 



' 4 (sin6 , e' fc m )$ fc , m + -^-e fc!m $' fe '. I 

sm 6 1 



Then ( p3| ) becomes 



Em=o^- m ^^ m ^ l .'» in if, 



(24) 



*U (°) = **.m («) = 



If we require that for all m 



-$km = Mm^fc.m in [0,0;] , 
*k, m (0) = Km («) = °, 



we obtain that $fc. m (y) = afc iTO cos (— if) satisfies (25) with fi r 
in (G4|) we have 



(25) 

1 2 

-tt- . Replacing this expression 



2_2 

I -r- 4 — 

fe,m i "t - sin 2 a 



1 7r m" r\ 



**, m (0) = $ 'fe,m («) = 0- 



Since the $fc jrn are independent, we have to solve, for every m, the Sturm-Liouville equation 



1 d 



{sm6Q' k<m ) + 



sin 9 d9 

Let us rewrite (^6|) in the following form 

1 d 



1 7r m 
sin 6 1 a" 1 



e fc , m = 0. 



fsinflei 



1 7r 2 m z 



sin 9 a z 



sin 9 d9 

so that we have to determine the eigenvalues \ k>m and the eigenfunctions of the operator 



(26) 



(27) 



1 d 



(sin 06' (0)) 



2™2 



1 ir z m 
sin 2 9 a 2 



-6(0). 



sin 9 d9 

In order to do this, let us consider the case a = n, that is the following equation 



1 cl 
sin 6d6 



[sm.6Q' k m ) + ; \ m 2 O fc , m = X k , m &k,m- 
sm 9 



(28) 



Now, for every m, ( pq ) has solution if Afe. m = fc (fc + 1), with fc > |m |, and the solutions are the 
Legendre polynomials @ kym (9) — P kym (cos9) (see for instance [Gro , [Ho], [Mu], [Mul]). Then, 
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for a given value of k, there are 2k + 1 independent solutions of the form Q k ,m (9) &k,m (<£>), one 
for each integer to with — k < m < k. Now, by the classical comparison principle, if we decrease a 
the corresponding eigenvalues X k . mi gi ven by (p7|), should increase, whereas if we increase a they 
should increase (see for instance [Cha|). More precisely, if to = the equations ( p7| ) and ( p8| ) are 
the same, therefore the eigenvalues do not change (and they are 0,2,6,...). If to > 1 we cannot 
give an explicit expression for the Xk, m for general a, but we can use the comparison principle. In 
conclusion, we obtain that each Y k = J2m=o ®k,m®k,m satisfies 

-A s *Y k = A fc , m F fc . (29) 

Now, one has that 

A (v k Y k ) = A r (v k ) Y k + -^v k A S 2Y k , (30) 

where _A_ r denotes_the Laplace operator in radial coordinates, that is A r = + Then, using 



fl29|) and (|30]), the condition (|2l|) becomes 

Afc ■, 



fc=0 



v k H — Ufc - pU p 1 v k 



Y k = 0. 



Since the Y k are independent, we get the following equations for v k : 

A Km (v k ) := -v'l - \' k + v k + ^v k - pU^Vk = 0, to = 0, 1, 2..., fc > 

Let us first consider the case to = 0. If k = 0, we have to find a vq such that 

2 

A),o («o) = -< - - w o + u o - pU^vq = 0. 



It has been shown in [Kw|, Lemma 6, that all the solutions of A 0j o (v) = are unbounded. Since 
we are looking for solutions Vq € W 1 ' 2 (R), it follows that vq = 0. 
For k = 1 we have to solve 

2 2 

-Al.O (Vi) = -v'l - -v[ +V! + -2«i - pU p ~ 1 vi = 0. 

Let (7 (r) denote the function such that U (x) = U (\x\), where U (x) is the solution of (p0|). 
Reasoning as in the proof of Lemma 4.1 in [AM], we obtain that the family of solutions of A\ q (ui) = 
0, with vi G W 1 ' 2 (M), is given by «i (r) — cU' (r), for some cel. 

Now, let us show that the equation A ky o (v k ) = has only the trivial solution in W 1,2 (R), 
provided that k > 2. First of all, note that the operator Ai o has the solution C7' which does not 
change sign in (0, oo) and therefore is a non-negative operator. In fact, if a denotes its smallest 
eigenvalue, any corresponding eigenfunction tp a does not change sign. If a < 0, then ip a should be 
orthogonal to U' and this is a contradiction. Thus a > and A\ t o is non-negative. Now, we can 
write 

A a _i_ A M ~ 2 

^fc,o — Ai.,0 H ^2 • 

Since A^^o — 2 > whenever k > 2, it follows that A^ q is a positive operator. Thus A k Q (v k ) = 
implies that v k = 0. 

If to > 1 and a < 7r, using the comparison principle, we obtain that each X ktm is greater than 
2. Then, reasoning as above, we have that each v k = 0. 
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Let us consider the case a > n. If m = 1 and k = 1, using again the comparison principle, 
we have that < Ai,i < 2; whereas for m = 1, k > 2, and for m > 2, fc > m, we have that 
each Afe. m > 2. Then in the last two cases we can use the non-negativity of the operator and 
conclude that Vk — 0. In the case m — 1 and k = 1 we note that the operator 

Ax,! := -v'{ - \\ + V! + ^±vr - pUP- 1 v 1 

has a non-negative eigenvalue, instead of the eigenvalue 0, since Ai.i < 2. Then also v\ = 0. 

Putting together all the previous information, we deduce that any v G Ker [/" (U)] has to be of 
the form 

v (x 1 ,x 2 ,x 3 ) = cU' (r) Yi (9, ip) . 

Now, Y\ is such that —Ag^Yx = Xi im Yi, namely it belongs to the kernel of the operator — A52 — 
Ai, m W, and such a kernel is 1-dimensional. In conclusion, we find that 

v G span ^U'Yi^ = span /— — 1 = TjjZ. 

This proves that i) holds. It is also easy to check that the operator 1'^ (U) is a compact perturbation 



of the identity, showing that ii) holds true, too. This complete the proof of Lemma 3.1. □ 

Remark 3.2. Since U is a Mountain-Pass solution of (20), the spectrum of 1% (U) has one 
negative simple eigenvalue, 1 — p, with eigenspace spanned by U itself. Moreover, we have shown 
in the preceding lemma that A = is an eigenvalue with multiplicity 1 and eigenspace spanned by 
J^-. If a < 7r the rest of the spectrum is positive. Whereas if a > n there is an other negative 

simple eigenvalue, corresponding to an eigenfunction U given by 



U (r, e,Lp)^u (r) cos (-<p) 6 (<9) . 



where & satisfies (|26|) with m = 1 and k = 1, and ii satisfies the equation 

-v" - ?-v' + v + ^-v-pU p - 1 v = 0. (31) 

From ( pl| ) one has that there exists a positive constant C such that, for r sufficiently large, u (r) < 
Ce~ r I. In conclusion, one has the following result: 



Corollary 3.3. Let U and U be as above and consider the functional Ik given in (p_9[). Then for 
every x\ G M, U Xl (x) — U (x — [x\, 0, 0)) is a critical point of Ik- Moreover, the kernel of I K (U) 
is generated by . If a < it the operator has only one negative eigenvalue, and therefore there 
exists S > such that 

dU 

I K (U)[v,v] ><S|M| 2 , for every v £ W 1 ' 2 (K) ,v±U, —. 

If a > it the operator has two negative eigenvalues, and therefore there exists 6 > such that 

dU 

I' K (U) [v, v]>5 \\v\\ 2 , for every v G W ia (K) , v±U, U, — . 
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4 Proof of Theorem 

For every Q on the edge T of <9f2 £ , let fi — min {/Xj}, so that in Bit (Q) f] il e we can use the new 
set of coordinates z. Now we choose a cut-off function (p^ with the following properties 

<Pn (x) = 1 in Bit (Q) , 

<p„ (x) = in R 3 \ Bsl (Q) , (32) 
|V<^| + |VVI<C in Bit (Q)\ Bit (Q) , 

and, for any Q £ T, we define the following function, in the coordinates (zi, z 2 , Z3), 

L/ Q , e (z) = <^(ez)[/(z-Q), (33) 

and the manifold 

Z e = {C/ Q , e : Q e T} . 

Now, we estimate the gradient of I e at Uq^ 7 showing that Z e constitute a manifold of pseudo-critical 
points of I € . 

Lemma 4.1. There exists C > such that for e small there holds 

\K(U Q ,e)\\ <Ce, forallQer. 



Proof. Let v <E W 1 ' 2 (^ e )- Since the function Vq^ is supported in B := Bj± (Q), see (p3[), we can 
use the coordinate z in this set, and we obtain 

I'e(U Q ,e)[v} = [ ^±vda + f (-A- g U Q , e + U Q , £ -\U Q>e nvdV s (z) 
= I + 11. 

Let us now estimate /: 



I = / - vda 1 + / — g vda 2 = 1 1 + h- 
Jdn el ov x J dQt2 dv 2 

If K = K a (Q) denotes the cone of angle equal to the angle of the edge in Q, we have 



OK 



h = / (Uq (z) Vip,i (ez) -vi + ipii (ez) VUq (z) ■ i7i )vddi 

JdK 

U Q (z) (ez) ■ (e (A Q (z l ,z 2 )) + e 2 D Q {z u z 2 ) , -1 + \<?\M z 2 ) | 2 

+ip f , (ez) VU Q (z) ■ (A Q (z l5 z 2 )) + e 2 D Q (z u z 2 ) , -1 + \^\ A Q (^2) | 2 

v(l + 0(e 2 |(z 1 ,z 2 )| 2 ))dz 1 dz 2 
= a + b. 

Since (e-) is supported in M 3 \ B^_ (Q) and Uq has an exponential decay, we have that, for e 
small, 

\a\ < Cee~£ [ \v\dz x dz 2 . (34) 

J OK 
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On the other hand 
b = 



&<\z-Q\<i 



-Q\<i 



(Pa (ez) VU Q (z) ■ (e (Aq {z 1 ,z 2 )) + e 2 D Q (z u z 2 ) , -1 + \e 2 \A Q z 2 ) \ 2 ^j 

v (l + O (e 2 | (zi,z 2 ) | 2 )) d Zl dz 2 
<p M (ez) VUq (z) ■ (e {A Q ( Zl> z 2 )) + e 2 D Q ( Zl , z 2 ) , -1 + \^ 2 \A Q z 2 ) \ 2 ^j 

v(l + 0(e 2 \{ yi ,y 2 )\ 2 ))d yi dy 2 



<Cee~^ \v\d Zl dz 2 + Ce \VUq\ ■ \v\d Zl dz 2 . (35) 

JdK JdK 



The estimates (|3j) and (|35|), and the trace Sobolev inequalities imply 

|/i|<Ce|M|. 
In the same way we can estimate I 2 , getting 

\I\<Ce\\v\\. 

Now let's evaluate II. Using dlj) one has 



(36) 



K 



2 7 " (0)* -V^ 



H= [ (-AU Qtt + U Q , e -\U Qte \P)vdVg(z) 

J K 

dU Q , 



dzi 



+ h"(o)-v (Z2 , Z3) u Q , £ ) 



VdVg (z) 

+0(e 2 ) f (^^UQ.A + W^UQ^vdV^z) 

JK 

= Ih + elh + O (e 2 ) Ih. 

Since AUq^ = UqAlp^ (ez) + 2VC/q • V</? M (ez) + ip^ (ez) AUq and both Aip^ (e-) and Vy^ (e-) are 
supported in K 3 \ Bj±_ (Q) : we get 



■<\x-Q\<& 



■<|*-Q|<; 



(-UqAip^ (ez) - 2VU Q ■ V<^ (ez)) v (1 + O (e\z\)) dz 



(-^ (ez) AU Q + U Q , e - \U Q , e \ p ) v(l + (e\z\)) dz 



l*-QI<£ 



(-AU Q + Uq- \Uq\p) v(1 + (e\z\)) dz. 



(37) 



Since Uq is a solution in R 3 the last term in ( |37| ) vanishes, and using the exponential decay of Uq 
at infinity and the properties of the cut-off function (see ([32|)) one has 



\Ih \ < Ce~£ ( \v\dz 

JK 
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By (|33|) we can compute also V( Z2)Z3 ) and V (z 2 ,z 3 )Uq,e and we have 



JJa = / 2 7 " (0) z! 



V 



( Z 2,Z3) " 



U Q + V (22,23)</?M ( £Z ) 



-2 7 " (0) zi 



diff, (ez) 



(22,23) 



dUn 



£<|2-<3|<£ 



2 7 " (0) *i 



+7" (0) • [V (Z2:Z3) <^ (ez) C/ Q + ^ (ez) V (z2iZ3 )£/ Q ] udVg (z) 
+7" (0) • V (Z2 , Z3) ^ (ez) [/ Q «dV§ (z) 



z-Q\<i 



27" (0) *i • V (z2 ,, 3 )^f + 7" (0) • V {Z2 , Z3) U Q vdV s (z) . 



Hence 
\Ih\<C 



tL<\z-Q\<£ 



2h." i"i I • i-il ( \Uq\ + |^2| + |v ( ^ s) tf | j - b" c» ; • r,,| 



|«|dV& (z) 



k-QI<£ 



2|<^ (ez) I • sup | 7 " (0) | (\z^\ • |V ( , 2l<s) ^| + |V (z2 , Z3) (7q|) (z) . 



Using again the exponential decay of Uq at infinity one can estimate the first term by Ce ft |w|dz 
and conclude that the second term is bounded. In the same way we can estimate J/3, getting 



\II\<Ce\\v\ 
From (pq) and (B8h we obtain the conclusion. 



(38) 
□ 



Now, we need a result of non-degeneracy, which allows us to say that the operator I" (Uq. e ) is 
invertible on the orthogonal complement of Tjj q e Z e . 

Lemma 4.2. There exists S > such that for e small, if a < it, there holds 

or j 

l"(U Q , e ){v,v] >S\\vf, for every v G W 1 ' 2 (Q e ) ,v±U Q , e , Q ' e 



dQ ' 

dUg,, 
e dQ 



and, if a > it , f/iere /io/c?s 

I'J (U Q , e ) [v, v}>6 \\vf , for every v G VF 1 ' 2 (fi e ) , uL£/ , e , tfq. 

where Uq, e is defined as Uq, e in (|33|). 

Proof. Let us consider the case a < it. Let i? 3> 1 and consider a radial smooth function xr '■ 
M 3 ->■ R such that 



and we set 



Xr(x) = 1 inBfl(O), 
X fl (z)=0 in R 3 \ 5 2fl (0) , 
|Vxfl|<| in B 2R (0) \ B fl (0) 



"l (z) = Xfl (x-Q)v (x) ; u 2 (x) = (1 - XR (% - Q)) v (x) . 



(39) 



1G 



A straight computation yields 

H 2 = IN| 2 + M 2 + 2 / (Vvi ■ Vv 2 + Vl v 2 ) dx. 

We write J n (Vv± ■ Vt>2 + V1V2) dx — 71 + 72, where 

7i= / Xb.{1-Xr)(v 2 + \Vv\ 2 )dx; 

72 = / (v2^v ■ Vxr - viVv ■ Vx_r - v 2 \\7xb,\ 2 ) dx. 

Since the integrand in 72 is supported in B 2 r (Q) \ Br (Q), using ( p9[ ) and the Young's inequality 
we obtain that I72I = or (1) ||i>|| 2 . As a consequence we have 

\\v\\ 2 = \\v 1 \\ 2 + \\v 2 \\ 2 + 2 1i +or(1)\\v\\ 2 . 

Now let us evaluate I" (Uq^) [v, v) = o\ + a 2 + 03, where 

ai = I'J {U Qt6 ) [ Vl , Vl ] , a 2 = I'l (U Q , e ) [«a,«a] , a 3 = H'l {U Q>e ) [v 1} v 2 ] . 

Similarly to the previous estimates, since Uq decays exponentially away from Q, we get 

^ 2 >C'- 1 ||i; 2 || 2 + 0e , iJ (l)H 2 ; 

a 3 > CTS + (1) Ibf • (40) 

Hence it is sufficient to estimate the term o\ . From the exponential decay of Uq and the fact that 
v±Uq^, dgQ- g it follows that 

(Vl, ^Q,e) w i,2 ( n e ) = - i V 2, U Qje ) wli2(Qc) = O e>R (1) ||t>|| 2 , 

vJ-^A =-U*j9A =o e ,,(l)||,|| 2 . (41) 



Moreover, since Uq^ is supported in B := Bg_ (Q), see (|33j), we can use the coordinate z in this 
set, and we obtain 

( v i, u Q,e)w^ 2 m 1 = / v 1 dU ®' i vd5+ f ui (-A s U Q , e + U Qj(L ) dV s (z) 
'on, av Jn 



(v 1 ,U Q ) wl , 2{K) +o e (l)\\v 1 \\, (42) 



where K — K a is the cone of opening angle equal to the angle of T in Q. In the same way we can 
obtain that 

dU Q , e \ ( dU c 



= h'flO +«e(l)lki||. (43) 



From the estimates (|4lJ), (^2|) and fl43|), we deduce that for R sufficiently large and e sufficiently 
small 



(vi,Uq) w1 , 2{k) = o e , fl (l) ||ui|| , 
9Q 



t*,^l = "../, CD ' , 



17 



Now we can apply Lemma p.l| , getting 

l"(U Q )[v 1 ,v 1 ]>S\\v 1 \\ wl , HK) +o e , R (l). 
Then the following estimate holds 



"i = I"(U Q )[v 1 ,v 1 ]+o e (l)\\v 1 \\ 

> 8\\vi\\ w i,2( K ) +°£,-r( 1 ) IMI 

> *||«i||+o e , K (l)|H|. (44) 
In conclusion, from ( fj(i| ) and Q) we deduce 

i' £ '(u Q , £ )[v,v]>s\\v\\ + o £ , R (i)\\v\\> S -\\v\\, 

provided R is taken large and e sufficiently small. This concludes the proof. 

The case a > n has substantially the same proof, but we h ave to consider also the function U 
and use the exponential decay of u at infinity, see Remark |3.2| . □ 

The following lemma provides an expansion of the functional I e (£/q, £ ) with respect to Q. 
Lemma 4.3. For e small the following expansion holds 

h (U Q , e ) = C a (Q) + 0(e), (45) 

where 

^=G-^T)rr |c/Q(r)rvsin2 ^ 



Proof. Since the function Uq^ is supported in B :— (Q), see (33), we can use the coordinate 
z in this set, and we obtain 

I* W Q , e ) = \[ (\V~ g U Q A 2 + U%, € ) dV s {z)--±- f \U Q A P+l dV~ g (z) . 
L JBnn e ' P "+" 1 JBnn e 

Integrating by parts, we get 

z JBndn, ov 
+ ~ [ U Q , £ (-A s U Q , e + U Q , e ) dV s ( Z )--L- f \U Q , e \" +1 dV s (z) 

= 1 + 11. 

Now, / can be split in two terms which correspond to the surface integrals on the "faces" of the 
edge T: 

If dU Q , e If dU Q , £ 

z JBndn el av i z JBndn e2 av i 
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It is sufficient to evaluate I\, since the estimate of J 2 is similar. Using the expression of Uq^ (see 
@) we get 



h 



Bndn el 



U Q<e ([/ Q V</v (ez) + ip„ {ez) VU Q ) 



(A Q (z lt za)) + e 2 D Q {z x ,z 2 ) , -1 + \( 2 \Aq «i) | 2 ) (l + (e 2 | (z 1; z 2 ) | 2 )) dz x dz 2 

= \l (ez) E/qV^ (ez) 

* ./£<|*-Q|<£ 

^ (A Q (zi, z 2 )) + e 2 C Q (z l5 z 2 ) , -1 + \e 2 \A Q (z u z 2 ) \ 2 ^j (l + O (e 2 | (z 1; z 2 ) | 2 )) dz x dz 2 



\z-Q\<i 



ifl (ez) UqVUq 



■ (e (Aq (z 1; z 2 )) + € 2 Dq (zx, z 2 ) , -1 + ^e 2 \A Q (z x ,z 2 ) | 2 ^) (l + O (e 2 | (z 1; z 2 ) | 2 )) dz x dz 2 . 

Similarly to the previous estimates, we get 1\ = O (e~^) + O (e). Then we obtain that 

l=0{e). 

Now, we have to evaluate /J. 



11=1 
2 



Bnn e 



U Q , e (-AU Q>e + U Q , e ) (l + O (e|z|)) dz 



+ - 



Bnfi e 



Uq, 



27" (0) zt ■ V (z2 . Z3) -^ + 7" (0) • V (Z2 , Z3) U Q 



p+1 



Bnn, 



(l + 0(e\z\))dz 
+0(e 2 \ Zl \ 2 ) 
U Q , e \ p+1 (l + 0(e\z\))dz. 



We have 



II 



1 1 



2 J5+ 1 



a(Q) 



OO p7T 



\Uq (r) | p+1 r sin 2 fldrdfl + O (e) 



Putting together (ff6|) and (47), we obtain (Eft) and this concludes the proof. 



(46) 



(47) 
□ 



Let Pq : W 1 ' 2 (f2 e ) — ► (Ttf S( Z e ) be the projection onto the orthogonal complement of 
Tfj Q c Z e , for all Q on the edge T of dfl e . According to the lemma 4.2, we have that for e sufficiently 
small the operator Lq — Pq o I'J (E/q )£ ) o Pq is invertible and there exists C > such that 



L7} 



< C. 



Now, using the fact that I'J (Uq^) is invertible on the orthogonal complement of Tu Q t Z t , we will 
solve the auxiliary equation. 

Proposition 4.4. Let I t be the functional defined in (^J). Then for e > small there exists a 
unique w = w (e, Q) 6 (Tu Qe Z e ) such that I' e (Uq jC + w (e, Q)) € Tjj q e Z e . Moreover the function 
w (e, Q) is of class C 1 with respect to Q and there holds 

dw (e, Q) 



\\w(e,Q)\\ <Ce, 



dQ 



< Ce. 



(48) 
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Proof. We want to find a solution w G {T Uq € Z^ of PqI^Uq^ + w) — 0. For every w G 
(Tu Q € Z e ) we can write 

l[ (U Q , t +W)=I[ (U Q , e ) + I'J (U Q , e ) [W] + i?Q, e (W) , 

where Rq^ (w) is given by 

R Q , e (w) = I' e (U Q , £ +w)-I' e {U Q , e ) - I': (U Q , e ) [w] . 
Given v G W 1 ' 2 (Q e ) there holds 

%WH = - / (\U Q ,e+w\ P - \U Q , e \ P - P \U Q . e \ p - 1 w)vdx. 

Using the following inequality 

| (a + b) p -a? - V a p - X b\ < 
(a, b G R, \a\ < 1), the Holder's inequality and the Sobolev embeddings we obtain 

\\R Q ,e( W )[v}\\ < C [ (\ W \ 2 + \w\ P ) \v\dx 

< c(\\w\\ 2 + \\w\F)\\v\\. (49) 

Similarly, from the inequality 

\{a + b l f -(a + b 2 f -pa p - 1 {b l -b 2 )\ 



C (p) \b\P for p < 2, 
C (p) (\b\ 2 + \b\P) forp>2, 



< 



C(p) (l&i I 33 " 1 + HP' 1 ) \h - b 2 \ for p < 2, 
C(p) (|6i| + N + l^iT 1 + |6 2 r 1 ) \h - b 2 \ for p > 2, 



(a, &i, 62 £ K, |a| < 1) we get 

\\R Q , e { Wl )[v]-RQ, t (w 2 )[v}\\ 
< C [ (\ Wl \ + \w 2 \ + \wi\P- 1 + \W2\P- 1 ) \ Wl - w 2 \ ■ \v\dx 

<c(Kn + iKii + Kir i +i| W2 ir i ) 

K-itf2||-IH|. (50) 

Now, by the invertibility of the operator Lq = Pq o I" (Uq^) o Pq, we have that the function w 
solves PqI'c (Uq, e + w) = if and only if 

w = - (Lq)- 1 [P Q I> (U Q>e ) + PqRq.c (w)} . 

Setting 

N Q , e (w) = - (Lq)- 1 [P Q l' e (U Q , e ) + P Q R Q , e H] , 

we have to solve 

w = N Q ^ (w) . 
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The norm of I' t (E/q, c ) has been estimated in Lemma 1.1. Then from (^9|) and ( p0| ) we obtain the 
two relations 

\\N Q , e (w)\\<G\e + C 2 (\\w\\ 2 + \\w\\ p ); (51) 
\\N Qti (wx)-N Q , e (w 2 )\\ < 
C (\\w x \\ + \\w 2 \\ + \\wtf- 1 + \\w 2 \f- x ) |K — wall - (52) 

Now, for C > 0, we define the set 

W e = [w € (Tu^Ze) 1 - : \\w\\<Ce]. 

We show that Nq^ £ is a contraction in Wq for C sufficiently large and for e small. Clearly, by ([pl|), 
if (7 > 2Ci the set Wq is mapped into itself if e is sufficiently small. Then, if wi,iU2 £ VF^, by 
(|52] ) there holds 

||iV Q>e {wi) - N Q , e (w 2 )\\ < 2C (Ce + C^e^ 1 ) \\ Wl - w 2 \\ . 

Therefore, again if e is sufficiently small, the coefficient of \\wi — w 2 \\ in the last formula is less 
than 1. Hence the Contraction Mapping Theorem applies, yielding the existence of a solution w 
satisfying the condition 

IMI < C e - ( 53 ) 

This concludes the proof of the existence part. 

Now the C 1 -dependence of the function w on Q follows from the Implicit Function Theorem 
(see also [AM], Proposition 8.7). In order to prove the second estimate in (fl8|), let us consider the 
map H : R 3 x W 1 - 2 (fi e ) xlxl — > W 1 ' 2 (fi £ ) x E defined by 



H (Q,w,a,e) 



Then w G (Tu Qt€ Z e ) is a solution of PqI' € {Uq^ + w) = if and only if H (Q,w,a,e) 
Moreover, for v e W 1 ' 2 (O e ) and /3 G R, there holds 

dH _ sr ^ / ^'([/ Q , e + W )M-/3 ac ■■■ ' 



(Q,w,a,e) [u,/3] = I c * "'"/ du , ' dQ (54) 



d(w,a) ' ' ^ ^u, — 

/ / e '([7 Q>e ) M -/3^ 
= l (^) 



o(hi +iHr x ) 



To prove the last estimate it is sufficient to use the following inequality 



l[a + b) _a 1 " \ C(p)(\b\ + \b\P- 1 ) forp>2, 



(a,b G M, |a| < 1), the Holder's inequality and the Sobolev embedding. Using the invertibility 

dH 
d(w.a) 



of the operator Lq — Pq o I" (Uq^) o Pq, it is easy to check that d P^ a) (Q,0,0, e) is uniformly 



invertible in Q for e small. Hence, by ( |53| ) and (^4|), also Qp^ a \ {Q,w,a,e) is uniformly invertible 



in Q for e small. As a consequence, by the Implicit Function Theorem, the map Q 1— > (wq,C(q) is 
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of class C 1 . Now we are in position to provide the norm estimate of dW gq^ ■ Differentiating the 



equation 



with respect to Q, we obtain 



H (Q,w Q ,a Q ,e) = 



OH 
0= — 



dH d(w Q ,a Q ) 



d (w, a 



dQ 



Hence, by the uniform invertibility of (Q, w, a, e) it follows that 



d(wq,aq) 



< C 



dQ 

rj (u Qtt + w 

: C 



< C 
dUn 



dU Q „ 



dQ 

d 2 u Q , 



I' e '(U Q , e + w) 
d 2 U Q , 



dQ 

l' e ' (U Q , 6 + w) 



dQ 2 
dUn. 



dQ 



dQ 2 
+ \a\ + || 1X7 ]1 + e 



Note that a, similarly to w, satisfies \a\ < Ce. By the estimate in (M) we obtain 



dU, 



Q-< 



dQ 

OUq 
dUr 



i': (u Qtt + w) 

Using the fact that I" (Uq) 
I'J {Uq, e + w) 

For any v 6 W 1 ' 2 (K), one finds 

\(I' e '(U QiE )-I"(U Q )) 



< 



dQ 



i':(u Q , e ) 

we obtain 
I"{Uq 



dU, 



dQ 



+ c(|| U ;|| + || W |r 1 ) 



dUr 



dzi 



I"(Uq) 



dUq 
dzi 



+Ce + c(|| W || + ||«;||^ 1 ) 



dUq 
dz\ 



< 



p[ \U Q , e -U Q \^v + Ce. 
Jxnn r dz x 



The last three formulas implies the estimate for dW gQ^ ■ This concludes the proof. 



□ 



Now we can state the following result, which allows us to perform a finite-dimensional reduction 
of problem (^) on the manifold Z t . 

Proposition 4.5. The junctional ^ € : Z t — > R defined by ^ e (Q) = I e (Uq ye + w (e, Q)) is of class 
C 1 in Q and satisfies 



K (Q) = o 



r e (U Q , £ +w(e,Q)) = 0. 



Proof. This proposition can be proved using the arguments of Theorem 2.12 of AM]. From a 
geometric point of view, we consider the manifold 



Z e = {U Q , e +w{e,Q):Q€T}. 
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Since ( Eg ) holds, we have that for e small 

Tu Q ,,Z e ^ T UQ e+w ^ Q )Z e . (55) 

If Uq^+w (e, Q) is a critical point of I e constrained on Z e , then I' e {Uq i£ + w (e, Q)) is perpendicular 
to T(j Q e +w(e,Q)Zei an d hence, from (p5|), is almost perpendicular to T\j Q c Z t . Since, by construction 
of Z t , it is I' e {U Q . t + w (e, Q)) G T Uq e Z e , it must be J £ (U Q , e + w (e, Q)) = 0. This concludes the 
proof. □ 



4.1 Proof of Theorem [O 

First of all we have 



* e (Q) = I e (U Qt£ + w(e,Q)) 

= h (U QtE ) + I' e (U Q , e ) [w (e, Q)] +0(\\w (e, Q)\\' 



Now, using Lemma 4.1 and the estimate (^8|) we infer 

* e (Q) = / 6 (f/ Q ,e) + O (e 2 ) . 

Hence Lemma |4.3| yields 

f £ (Q) = C a (Q) + O (e) . 

Therefore, if Q S F is a local strict maximum or minimum of the function a, the thesis follows 
from Proposition |4.5| . 
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